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Abstract: 

We present a determination of the parton distributions of the nucleon from a global 
set of hard scattering data using the NNPDF methodology including heavy quark mass 
effects: NNPDF2.1. In comparison to the previous NNPDF2.0 parton determination, the 
dataset is enlarged to include deep-inelastic charm structure function data. We implement 
the FONLL-A general-mass scheme in the FastKernel framework and assess its accuracy 
by comparison to the Les Houches heavy quark benchmarks. We discuss the impact 
on parton distributions of the treatment of the heavy quark masses, and we provide a 
determination of the uncertainty in the parton distributions due to uncertainty in the 
masses. We assess the impact of these uncertainties on LHC observables by providing 
parton sets with different values of the charm and bottom quark masses. Finally, we 
construct and discuss parton sets with a fixed number of flavours. 
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1 Introduction 



The inclusion of effects related to heavy quark masses in the determination of parton 
distributions (PDFs) has received an increasing amount of attention over the last few 
years, driven by the increase in accuracy and reliability in the determination of the PDFs 
required for phenomenology at the LHC (see e.g. Ref. [Ij and references therein). Only a 
few years ago, PDFs in common use (such as e.g. CTEQ6.1 [2j) were based on the so-called 
zero- mass variable- flavor number scheme (ZM-VFN), in which heavy quarks decouple at 
scales below their mass, < but are otherwise treated as massless partons, which 
amounts to neglecting all contributions of order m\/Q'^. While this approximation only 
applies to heavy quark distributions in the vicinity of their respective thresholds, the 
ensuing modification of the initial conditions to perturbative evolution for the heavy quark 
distributions also affects light quark PDFs (the momentum sum rule means that a change 
in any quark's momentum fraction must be accompanied by corresponding changes in 
the momentum fractions carried by all other partons). The high-energy behaviour of 
light quark PDFs may then be affected at the level of several percent, and the ensuing 
shift in predictions for precise high-energy standard candles such as the W and Z cross- 
sections may be quite significant [3]. Furthermore, observables which depend on the heavy 
quark distributions (such as the single-top production cross-section, which probes the h 
distribution) are substantially affected 

The purpose of this paper is threefold: first, to present a determination of parton 
distributions based on the NNPDF methodology [51-110] with heavy quark mass effects 
included. Second, to provide tools to study uncertainties related to heavy quark masses 
and more general heavy quark effects. Third, to assess the impact of these uncertainties 
on phenomenology. 

The first goal will be achieved by repeating a next-to-leading order global PDF determi- 
nation based on exactly the same methodology used in the construction of the NNPDF2.0 
PDF set [10], but now with heavy quark mass effects included up to order as through 
the so-called FONLL-A scheme [llj. The FONLL method (first suggested in Ref. [S] 
and generalized to deep-inelastic scattering in Ref. [11]) is especially convenient in that it 
allows the inclusion of heavy quark mass effects to any desired order in and any desired 
logarithmic order. The FONLL-A version corresponds to the combination of 0(qs) mass 
effects with NLO evolution equations and coefficient functions: at this NLO-O(as) order, 
the FONLL method coincides with the so-cahed S-ACOT (simplified [13] ACOT [Ti] ) 
method, adopted for instance in the CTEQ6.6 [15] and CTIO [16j NLO PDF determina- 
tions. The dataset used here also coincides with that of NNPDF2.0, but supplemented by 
charm deep-inelastic -F| structure function data. 

The second goal will be achieved by providing sets of parton distributions which corre- 
spond to different values of the heavy quark masses: the uncertainty related to the choice 
of the quark mass can then be determined simply by variation of the mass value, while 
combined PDF+m/j uncertainties can be determined by constructing Monte Carlo sets 
of replicas in which the mass is varied according to a probability distribution (typically 
gaussian) with a suitable width, in analogy to what was done in Refs. |17pi8] to determine 
combined PDF+a^ uncertainties. We also provide a determination of the correlation be- 
tween the heavy quark masses and individual PDFs. Finally, we will provide PDF sets 
with various fixed number of flavours. 
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The third goal wih be achieved by computing PDF and heavy quark mass uncertainties 
for various LHC standard candles: W, Z, Higgs and top production. We will also present 
a preliminary estimate of theoretical uncertainties related to higher order heavy quark 
mass corrections. 

The outline of this paper is the following: in Sect. [2] we discuss the features of the 
datasets included in the NNPDF2.1 analysis, with emphasis on the ZEUS and HI data on 
the charm structure function. Then in Sect. [3]we review the FONLL scheme of Ref. pT] for 
the inclusion of heavy quark mass effects in neutral current structure functions and present 
its generalization to charged current deep-inelastic scattering. In Sect. |l]we present the 
NNPDF2.1 PDF set and compare it with previous NNPDF releases and with the other 
global PDF sets, while in Sect. [5] we perform the same comparisons for LHC standard 
candles, thus elucidating the impact of the inclusion of heavy quark mass effects in the 
NNPDF framework. In Sect. Owe explore the impact of the uncertainty on the values of the 
heavy quark masses both on PDFs themselves and on LHC processes using NNPDF2.1 sets 
with varying m/j. Finally, in Sect. [7] we present NNPDF2.1 sets with various fixed number 
of flavours. Technical details on the implementation and benchmarking of FONLL neutral 
and charged current structure functions in the FastKernel computational framework of 
Ref. [IQ] are collected in Appendices [A] and [Bj 
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2 Experimental data 



In this Section we discuss the experimental data used for the NNPDF2.1 analysis. First of 
all we motivate the kinematical cuts that are applied to our dataset. Then we present the 
details and kinematical coverage of the NNPDF2.1 dataset, with special emphasis on the 
new charm structure function data. Finally, we discuss the implementation of positivity 
constraints. These data have been used to generate Monte Carlo replicas, which have been 
checked to reproduce the statistical features of the original dataset. The replica generation 
and its testing has been performed in the same way as in previous NNPDF analyses [HlllOj 
and will not be discussed further here. 

2.1 Kinematical cuts 

The NNPDF2.1 dataset has been subjected to some kinematic cuts: specifically, the cut 
in W'^ is the same as in previous NNPDF fits, VF^in = 1^.5 GeV^, but the cut in is 
slightly higher. While in NNPDF2.0 the cut in for the DIS data was set to be Q^^^ = 2 
GeV^, in the NNPDF2.1 analysis we use a somewhat more restrictive kinematical cut in 
Q^, namely Q^i^ = 3 GeV^. There are two main motivations for this modification which 
we now discuss. First, very close to the heavy quark threshold the predictions for F| from 
the GM scheme might suffer from instabilities due to the threshold behaviour. One would 
like to avoid having data crossing the charm mass threshold when varying the heavy quark 
mass in various fits. This suggest to use a value of Q^j^ at least as large as the maximum 
value of the charm mass than can be considered acceptable. Q^jn = 3 GeV^ is then a 
reasonable choice since then m™*" ~ 1.7 GeV. Furthermore, there is now an indication of 
possible deviations from NLO DGLAP in the small-x and HERA data [HllO]. These 
deviations are mostly relevant in the smaller bins of HERA data. The theoretical 
uncertainty in the PDFs and LHC observables related to their inclusion in the global fit 
is moderate as compared to the PDF errors and other uncertainties, but removing the 
HERA points below Q^i^ reduces these theoretical uncertainties even further. The price 
to pay for this reduced theoretical uncertainty is an increase in statistical uncertainty: 
indeed, we will see in Sect. 14.31 that removing the data below Q^j^ = 3 GeV^ results in an 
increase of PDF uncertainty in the small-x gluon PDF due to the reduced experimental 
information. 

On top of the previous general kinematical cuts, applied to all the DIS experiments, 
we will also perform additional cuts on the HERA F2 data. The motivation for these is 
that in this work we will use the FONLL-A general-mass scheme for heavy quarks, and 
as discussed in , FONLL- A0 provides a poor description of the data in the smallest x 
and bins due to missing large O (a^) corrections. Only the FONLL-B scheme can cure 
this problem since it includes consistently O (a^) corrections in F2 into a NLO fit, as can 
be seen in [11] and we will review in Sect. [3j We will thus remove from the fit HERA F2 
data with < 4 GeV^ and data with < 10 GeV^ for x < 10~^. These cuts ensure 
that all F2 experimental data included in the fit are well described by O (as) theory. 

In Table [T] we summarize the choices for the initial evolution scale and kinematical 
cuts applied in this work, compared to the choices in other recent PDF determinations. 

^Note that this is true for any heavy quark scheme that does not include the O (q^) corrections, hke 
for example the S-ACOT-x used in the CTEQ/CT family of PDF sets. 
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Ql [GeV^] 






NNPDF2.1 


2.0 


3.0 


12.5 


NNPDF2.0 ilOj 


2.0 


2.0 


12.5 


CTIO fT6] 


1.69 


4.0 


12.25 


MSTW08 [21J 


1 


2.0 


15.0 


ABKM09 [22] 


9 


2.5 


3.24 


HERAPDFl.O [23] 


1.9 


3.5 


155.75 



Table 1: The values of the mitial evolution scale where the PDFs are parametrized, Qq, and the 
kinematical cuts in and applied to the fitted DIS dataset, Q^^jj, and W,^;^, in the present 
work and in other recent PDF determinations. As discussed in the text, further cuts are applied 
to F2 data in the NNPDF2.1 case. For HERAPDF the value of W^i,j given is the minimum of the 
HERA dataset and no cut is performed. 

Note that HERAPDF does not perform a cut in VF^ since they only include HERA data 
which do not extend to the low W"^ region. 

2.2 NNPDF2.1 dataset 

Now we discuss the datasets that are included in the present analysis. As compared to the 
NNPDF2.0 analysis [lOj . on top of all relevant data from DIS, Drell-Yan and weak vector 
boson productiorll and inclusive jet production we include here all the relevant charm 
structure function F2{x,Q'^) data from the HI and ZEUS experiments at HERA |25H31] . 
These datasets provide a handle on the small- a; gluon, and are sensitive also to the value 
of the charm mass rric- On the other hand, HERA F2 has much larger uncertainties, 
and is thus not included in the present analysis. The kinematical coverage of all the 
datasets included in NNPDF2.1 is summarized in Table [2] and in Fig. [TJ Note that the 
only differences with respect the NNPDF2.0 dataset are the addition of HERA F2 data 
and the new kinematical cut Q^i^ = 3 GeV^. 

Now we describe in turn the features of the various datasets included in 

the present analysis. For most experimental sets the full correlation is not available and 
thus one is forced to add in quadrature systematic and statistical uncertainties. The full 
correlation matrix for all data points, including the cross-correlations between datasets and 
between HI and ZEUS will be provided together with the combined HERA F2 dataset: this 
combination will thus significantly improve the accuracy of the existing separate datasets. 

The F2 data which we use in the NNPDF2.1 analysis are the following: 

• The ZEUS 96-97 D*^ analysis [25]. 

In this analysis F2 is extracted from the measurement of D*^ mesons reconstructed 
via their hadronic decays using data collected in the 1996 and 1997 running periods. 

• The ZEUS 98-00 D* analysis [26]. 

As in the previous case, F2 is extracted from the measurement of D*^ mesons 
reconstructed via their hadronic decays, and uses data collected in the running period 
between 1998 and 2000. 

^The impact of the leptonic W asymmetry data from the Tevatron, not included in NNPDF2.0, as been 
studied in Ref. [24] using the Bayesian reweighting technique. 
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Deep-Inelastic scattering 


— Experiment 


Set 


Ref. 


-'^dat 


^min 


^max 






NMC-pd 






260 (132) 












NMC-pd 


1321 


260 (132) 


0.0015 (0.008) 


0.68 


0.2 (3.5) 


99.0 


NMC 






288 (221) 












NMC 


1331 


288 (221) 


0.0035 (0.009) 


0.47 


0.8 (3.2) 


61.2 


SLAG 






422 (74) 












SLACp 


1341 


211 (37) 


0.07 (0.1) 


0.85 (0.55) 


0.58 (3.0) 


29.2 




SLACd 


1341 


211 (37) 


0.07 (0.1) 


0.85 (0.55) 


0.58 (3.2) 


29.1 


BCDMS 






605 (581) 












BCDMSp 


EH 


351 (333) 


0.07 


0.75 


7.5 


230.0 




BCDMSd 




254 (248) 


0.07 


0.75 


8.8 


230.0 


HERAI-AV 






741 (592) 












HERAl-NCop 


l23l 


528 (379) 


6.2 10" ' (4.3 10"'') 


0.65 


0.045 (3.5) 


30000 




HERAl-NCom 


nn 


145 


1.3 10"" 


0.65 


90.000 


30000 




HERAl-CCep 


[23] 


34 


0.008 


0.4 


300.0 


15000 




HERAl-CCem 


m 


34 


0.013 


0.4 


300.0 


30000 


CHORUS 






1214 (862) 












CHORUSnu 


mi 


607 (431) 


0.02 (0.045) 


0.65 


0.3 (3.0) 


95.2 




CHORUSnb 


l37l 


607 (431) 


0.02 (0.045) 


0.65 


0.3 (3.0) 


95.2 


FLH108 






8 












FLH108 


nn 


8 


0.00028 


0.0036 


12.0 




NTVDMN 






90 (79) 












NTVnuDMN 


I39II40I 


45 (41) 


0.027 


0.36 


1.1 (3.1) 


116.5 




NTVnbDMN 




45 (38) 


0.021 


0.25 


0.8 (3.1) 


68.3 


ZEUS-H2 






127 












Z06NC 


41 


90 


5 10"-* 


0.65 


200 


3 lO'' 




Z06CC 


42 


37 


0.015 


0.65 


280 


3 10=" 








HERA charm structure function 


data 










ZEUSF2C 






69 (50) 












ZEUSF2C99 


mi 


21 (14) 


5 10"=" (SIO"*) 


0.02 


1.8 (7.0) 


130 




ZEUSF2C03 




31 (21) 


310"'" (1.810"^) 


0.03 


2.0 (7.0) 


500 




ZEUSF2C08 




9 (7) 


2.2 10"* ( 6.5 10"*) 


0.032 


7.0 


112 




ZEUSF2C09 


nn 


8 


8 10~* 


0.03 


30 


1000 


H1F2C 






47 (38) 












H1F2C01 




12 (6) 


5 10"" 


3.2 10"" 


1.5 (12) 


60 




H1F2C09 




6 


2.4 10"" 


0.025 


120 


400 




H1F2C10 


mi 


26 


2 10"* (3.2 10"*) 


0.05 


5.0 (12) 


2000 


Fixed Target Drell-Yan production 


Experiment 


Set 


Ref. 


^Vdat 


f/^inin' a/^max 


[^min ■ ^max] 




A^max [GeV^] 


DYE605 






119 












DYE605 


m 


119 


1-0.20, 0.40J 


[0.14, 0.65J 


50.5 


286 


DYE866 






390 












DYE866P 


|44"45' 


184 


|0.0, 0.78J 




0.017, 0.87 




19.8 


251.2 




DYE866r 


46 


15 


[0.05, 0.53] 




0.025, 0.56 




21.2 


166.4 


Collider vector boson production 


Experiment 


Set 


Ref. 


^Vdat 


[ymin, !/max] 


[^min , ^max] 


A/,i„, [GeV-] 


M4,„. [GeV-] 


CDFWASY 






13 












CDFWASY 


[47l 


13 


[0.10, 2.63] 


[2.9 10"-', 0.56 


6463 


6463 


CDFZRAP 






29 












CDFZRAP 


l48l 


29 


[0.05, 2.85] 




2.9 10"", 0.80 




8315 


8315 


DOZRAP 






28 












DOZRAP 


[49] 


28 


[0.05, 2.75] 




2.9 10"", 0.72 




8315 


8315 


Collider inclusive jet production 


Experiment 


Set 


Ref. 


^Vdat 


[Hmin, !/max] 


[^min • J^max] 


PT.min [GcV^l 


Pr.max IGcV-] 


CDFR2KT 






76 












CDFR2KT 


50 


76 


[0.05, 1.85] 


4.6 10"", 0.90 


3364 


3.7 10^ 


D0R2CON 






110 












D0R2CON 


I5T1 


110 


[0.20, 2.20] 


|3.1 10"", 0.97 


3000 


3.4 10^ 


Total 


Experiment 






^Vdat 




Xmax 


Qi„„ [Gov-] 


<5max [Gov-] 


TOTAL 






4520 (3338) 


3.1 10"^ 


0.97 


2.0 


3.7 10-' 



Table 2: Experimental datasets included in the NNPDF2.1 global analysis. For DIS experiments 
we provide in each case the number of data points and the ranges of the kinematical variables 
before and after (in parenthesis) kinematical cuts. For hadronic data we show the ranges of parton 
X covered for each set determined using leading order parton kinematics. Note that hadronic data 
are unaffected by kinematic cuts. The values of Xmin and Q^jin for the total dataset hold after 
imposing kinematic cuts. 
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NNPDF2.1 dataset I 



10^ 



010" 

10k 

1 

10- 





X 


NMC-pd 






NMC 




+ 


SLAC 




+ 


BCDMS 






HERAI-AV 






CHORUS 




• 


FLH108 




• 


NTVDMN 




X 


ZEUS-H2 






ZEUSF2C 




O 


H1F2C 






DYE605 




+ 


DYE886 






CDFWASY 






CDFZRAP 






DOZRAP 




X 


CDFR2KT 






D0R2CON 



■ ^ w \^ 

y/cxK^''^ ' tIi ?R 




10 



10" 



10 



r2 



10-^ 



Figure 1: Experimental datasets which enter the NNPDF2.1 analysis. The kinematical coverage 
of each dataset is summarized in Table [H 



• The 04-05 ZEUS D^,D^ analysis [27]. 

In this analysis, based on the HERA-II running period of 2004 and 2005, D mesons 
are reconstructed via their hadronic decays. An improved precision is obtained 
reducing the combinatorial background to the D meson signals by using the ZEUS 
micro-vertex detector to reconstruct displaced secondary vertices 

• The 2005 ZEUS muon analysis [28]. 

This dataset is based on the measurement of muons that are generated in charm 
production from their semileptonic decays. Data was collected during the 2005 
HERA-II running period. 

• The HI 96-97 D*^ analysis |29]. 

This analysis, based on the 1996-1997 running period, used similar reconstruction 
strategies as the corresponding ZEUS analysis, namely the reconstruction of D*^ — )■ 
D^Tr~^ using the D* — mass difference method. 

• The HI large 04-07 D*^ analysis [30]. 

This analysis determines F2 via identified D mesons produced at large virtualities 
> 100 GeV^, and is based on data collected in the HERA-II running period 2004 
and 2007. 
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• The HI 1ow-Q2 06-07 D*^ analysis [HI] 

This is analogous to the previous measurement, but now covering the small and 
medium region. It is based on data obtained in the HERA-II 2006-2007 running 
period. Events containing heavy quarks are distinguished from those containing only 
light quarks using variables that are sensitive to the longer lifetimes of heavy flavour 
hadrons, like the transverse displacement of tracks from the primary vertex. 

There are more published F2 datasets from HERA but the ones that are included here 
supersede previous obsolete measurements and are the basis of the combined HERA F2 
dataset. In Sect. H]we will quantify the impact of the HERA F2 data onto the PDFs. 

A concern with F2 data which has been sometimes used to motivate their exclusion 
from PDF determinations is the fact that the way F2 is usually defined experimentally, 
as the contribution to F2 with at least one charmed quark in the final state, is affected 
by mass singularities (i.e., it is not finite in the limit in which rric — )• 0). Here we will 
adopt a definition of F2 (as the contribution to F2 when only the charm electric charge is 
nonzero) which is free of mass singularities; the deviation between this definition and that 
which is used to define the experimental observable is estimated in Ref. by means of 
a suitable resummation method, and shown to be negligible in the region of the HERA 
data. Also, F2 is affected by theoretical uncertainties related to the extrapolation from 
the experimentally accessible region (restricted in pT and rj) to the full phase space. This 
theoretical uncertainty is estimated using QCD exclusive partonic calculations and added 
as an extra source of systematic uncertainty in the experimental analysis. 

2.3 Positivity constraints 

As discussed in [10], within the NNPDF framework general theoretical constraints can be 
imposed guaranteeing that the fitting procedure only explores the subspace of acceptable 
physical solutions. An important theoretical constraint is the positivity of physical cross- 
sections. As discussed in Ref. [52], positivity should be imposed on observable hadronic 
cross-sections and not on partonic quantities, which do not necessarily satisfy this con- 
straint (except at leading order where the probabilistic interpretation holds). Positivity 
constraints may be implemented in various ways; here we will impose them through La- 
grange multipliers, i.e. in practice by adding pseudo-datasets for physical cross sections 
with extremely small uncertainties in such a way that negative cross sections would lead 
to a very large contribution to the x^- 

In NNPDF2.1 we impose positivity of the following observables: 

• The longitudinal structure function Fl{x,Q'^), which constrains the gluon positivity 
at small-x. 

• The charm production cross section in neutrino DIS, d'^a'^'^ / dxdy [9j, which con- 
strains the strange PDFs both at large and at small-x, beyond the reach of existing 
data. 

• The neutral current DIS charm structure function F2{x,Q'^), useful to impose the 
positivity of the gluon at very large-x, where it is not constrained by any experi- 
mental dataset. 
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All the positivity constraints are implemented at a low scale Qpos that we take to 
be Ql^g = 2 GeV^, in the range x G [lO"'^ 

) s^maxl 7 where Xmax is the corresponding 
kinematical boundary, Xmax ~ 0.1 for NC scattering and Xmax ~ 0.5 for CC scattering. 
DGLAP evolution then takes care or preserving the positivity properties for higher scales. 
We note that the physical observables for the pseudo-data that implement the positivity 
constraints are computed consistently at the same perturbative order as all other physical 
observables, in the present case next-to-leading order perturbative QCD. 
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3 Structure functions with heavy quark mass effects 



The FONLL-A general-mass scheme was introduced for neutral current structure functions 
in Ref. [11] . We begin this Section with a brief review of this scheme, emphasizing the 
impact of heavy quark effects on DIS structure functions. We then discuss the values of 
the heavy quark masses and the associated uncertainties adopted in the present analysis. 
The corresponding analysis for charged current structure functions is presented in the last 
part of this Section. 

3.1 The FONLL-A General Mass scheme for NC structure functions 

The FONLL general-mass scheme, originally proposed in the context of heavy quark 
photo- and hadro-production, was generalized in Ref. [TT] to deep-inelastic structure func- 
tions. We refer the reader to Ref. ^llj for a detailed discussion of the scheme, and for the 
notation adopted in this Section. The FONLL approach allows for a consistent combi- 
nation of terms determined in a massive, or decoupling, or fixed-flavour number (FFN) 
scheme, in which the heavy quark is subtracted at zero momentum (rather in the MS 
scheme), so it decouples for scales much below its mass, and it is included in Feynman 
diagrams up to some fixed order in above its threshold for its production, with terms 
determined in a massless, or zero-mass (ZM), or simply MS scheme, in which the heavy 
flavour is treated as another massless parton, so it is included in the all-order resumma- 
tion of collinear logarithms, up to a suitable chosen logarithmic order (LO, NLO, etc). A 
significant feature of the FONLL method is that the fixed perturbative order of the FFN 
computation and the resummed logarithmic order of the ZM computation which are being 
combined can be chosen independently of each others. 

In the present analysis, we combine FFN massive terms up to order as with a NLO ZM 
computation; this is called FONLL-A in Ref. [llj. As shown in Ref. [53], this turns out to 
be identical to the S-ACOT [13] scheme used in recent CTEQ/CT PDF determinations \15\ 
I16j . Once a specific "general mass" (GM) scheme for the combination of FFN and ZM 
terms has been chosen, there is still a freedom in the treatment of subleading terms: indeed, 
it turns out to be phenomenologically convenient to suppress subleading terms near the 
quark threshold (see Ref. |lHI54j). In this work we adopt the so-called threshold or 
damping factor method of Ref. ^T] for the treatment of subleading terms. In Ref. ^3] the 
damping factor method is benchmarked against various implementations of the alternative, 
commonly used x~scaling method for the treatment of subleading terms. 

We now present the explicit expressions for the F2^h heavy quark structure function H. 
The FONLL-A heavy quark structure function is given by the sum of two terms: 

fT'^H^, Q') = F^{x, Q') + e{Q'- ml) (l - ^) ' Q'y (1) 

The first contribution on the right-hand side of Eq. ([1]) is the massive-scheme heavy-quark 
structure function at O (a^): 

F^\x,Q') = X /' ^Cg) ^,asiQ')) 9^^^^'Kv,Q') • (2) 
Jx y \y '"'h / 

■^See Ref. [11] for the discussion on the FONLL expressions for the longitudinal structure functions. 
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The heavy quark gluon coefficient function is given by 



mt 

The O (os) coefficient is 



1 + ?; 



^'^''f^Al = e{W'-Aml)xTn[{z' + il-zf + Aezil-3z)-8e'z')log\- 

V "^h/ 1 — 

+{8z{l - z) - 1 - 4ez(l - ] , (4) 

where we have defined 



(5) 



and the partonic center of mass energy W"^ = Q^{\ — z)/ z. 

The second term on the right-hand side of Eq. ([1]) is the "difference" contribution 



X 



y 



a 



a 



2,g 

(ni + l) / X 



(ni+l) X 



y 



2,3 



y J 9''^ \y ml 



/i("'+l)(y,Q2) + /,K+l)(y^g2) 



+ 



(6) 



where h,h are the heavy quark parton distributions; at ffi'st-order in q^, B^^\ is given by 



9, h 



B 



(o),i 

9, 



^ <^^^ _ o„2^("i,0),l 



mt 



and the massless limit of the massive coefficient function is 



a 



(ni,0),l 
2,9 



z,— ]=Tr [z z) ) log 



mt 



m^z 



+ (8z(l -z)-l) 



(7) 



(8) 



which in the limit = m| reproduces as required the usual massless scheme coefficient 
function. 

Note that in all terms in Eq. ([1]) PDFs and are expressed in the same factorization 
scheme, namely, the decoupling Uf = 3 scheme. Exploiting this fact, it is easy to check 
explicitly that the "difference" term Eq. ([6]) is formally of higher order near the heavy 
quark threshold [11] (and thus in particular it can be suppressed using a suitable threshold 
prescription). 

Eq. ([1]) interpolates smoothly between the massive scheme at small and the massless 
scheme suitable at large Q^. As an illustration of the differences between various schemes 
for the heavy quark structure functions, in Fig. [2] we compare the i<2,c and the F^^c charm 
structure functions for various schemes: ZM, FONLL-A and the FFN scheme as a function 
of for different values of x. It is clear that FONLL-A interpolates smoothly between the 
FFN scheme near threshold and the massless scheme at large (also thanks to the use of 
a damping factor in Eq. ([TJ. For this comparison, PDFs and other settings, like the value 
of rric, are identical to those of the Les Houches heavy quark benchmark comparison 
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Figure 2: The charm structure functions F2^c{x,Q^) and Fl.c{x,Q^) as a function of for 
different values of x from x = 10~^ to x = 10~^ in various heavy quark schemes, computed using 
the FastKernel method: FONLL-A, ZM-VFN and the FFN scheme. The PDFs and settings are 
identical to those of the Les Houches heavy quark benchmark comparison. 
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ZM / GM Rel Diff in F2c ZM / GM Rel Diff in F2p 




10"^ 10""* 10"'' 10"^ 10"' 10"^ 10"^ 10"^ 10"^ 10"' 




10"^ lO"'* lO""* 10"^ 10"' 10"^ 10"'' 10"^ 10"- 10"' 



X X 

Figure 3: Left plots: the relative difference for the charm structure functions F2^c{x,Q^) and 
FL,cix,Q^) computed in the ZM and FONLL-A schemes as a function of a; and Q^. The PDFs 
and settings are identical to those of the Les Houches heavy quark benchmark comparison. Right 
plots: the same but now for the inclusive structure functions and F^{x,Q^). 

The comparison for the longitudinal structure function Fl^c shows that mass effects are 
much larger than in -F2,c; so the ZM computation is completely unreliable. 

The impact of heavy quark mass effects in DIS structure functions is further quantified 
in Fig.[3l where the relative difference between the ZM and FONLL-A schemes is computed 
as a function of x and Q^, both for the inclusive structure functions F2 and Ff^ and for 
the charm structure functions i<2,c and Fl^c- For the phenomenologically more relevant 
case of F2, we see that heavy quark mass effects can be as large as ~ 10%, decreasing fast 
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1+(F0NLL-B - FONLL-A) / FONLL-A 




10"^ lO-"^ 10-3 10"^ 

X 

Figure 4: Relative difference for the charm structure functions F2.c{x, Q^) between the FONLL-B 
and FONLL-A general- mass schemes, in units of FONLL-A, as a function of x and . The PDFs 
and settings are identical to those of the Les Houches heavy quark benchmark comparison. 

for increasing x and Q^. As in the case of Fig. [2] the Les Houches heavy quark benchmark 
settings have been used. Note that while the qualitative features of Fig. [3] are general, 
the quantitative detail can depend on specific features of the general-mass heavy quark 
scheme, like for example the prescription to suppress the subleading threshold terms. 

As discussed in Ref. it is possible to account for the phenomenologically relevant 
O (a^) corrections to i<2,c into an NLO PDF fit by means of the FONLL-B scheme. We 
show in Fig. mthe relative difference between i<2,c computed in the FONLL-B and FONLL- 
A schemes. As shown in the figure, their difference at small x and is rather large. The 
inadequacy of O (a^) theory to describe the low x and -F2,c data motivates the cuts to 
the -F2,c HERA datasets discussed in Sect. 12. ll At larger values of x and the differences 
between the two schemes become of the order of a few percent, much smaller than the 
typical experimental uncertainties, thus validating the inclusion of the i<2,c data into the 
present fit based on the FONLL-A scheme. 

The O {as) massive scheme heavy quark coefficient function, Eq. (I3]K was first com- 
puted in Refs. |55H57j . while its Mellin transform, hitherto not availabl^j is presented in 
Appendix |X1 Details of the implementation of the FONLL-A scheme in the FastKernel 
framework used in the NNPDF analysis are also given in Appendix [Al We have assessed 
the accuracy FONLL implementation computing the Les Houches heavy quark benchmark 
tables [53], showing that the accuracy is sufficient for precision PDF determination. 

numerical parametrization of the Mellin space heavy quark coefficient functions up to O (a^) was 
provided in Ref. [58] . 
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3.2 FONLL Charged Current structure functions 

The FONLL method for charged currents was only mentioned briefly in Ref. [11]: here 
we provide a detailed explanation of this general-mass scheme for CC structure functions. 
In the charged current sector, coefficient functions are only known up to 0{as), so the 
FONLL-A scheme can be constructed, while the FONLL-B and C which are based on 
unavailabl^ C(a^) massive results cannot. 

Heavy quark mass effects are required to describe charm production in neutrino DIS 
(the dimuon process) and to a lesser extent also the inclusive neutrino reduced cross 
sections, since in both cases most of the data lie close to the charm threshold, > m^. 
HERA charged current data on the other hand are at large and thus for practical 
purposes any general-mass scheme reduces to the ZM-VFN scheme. 

Here we generalize the FONLL-A scheme to charged current structure functions. Its 
implementation in the FastKernel framework requires the analytic computation of the 
Mellin transforms of the O (as) charged current heavy quark coefficient functions [61]. The 
detailed description of the implementation is given in Appendix [Bj We then benchmark 
the FONLL implementation in FastKernel, using the Les Houches heavy quark benchmark 
settings [53] comparing with a x-space code written for this purpose that implements 
FONLL-A for CC structure functions. For simplicity, we will make the assumption that 
\Vcs\ = 1, and the rest of the CKM matrix elements are zero. The generalization to realistic 
CKM elements, as actually implemented in FastKernel, is straightforward. We assume also 
a single heavy quark, the charm quark with mass rric- The factorization scale is set to 
be equal to /u|. = Q^. Finally, we consider only neutrino induced charm production, the 
anti- neutrino case is again straightforward. 

In the FFN massive scheme, the charged current charm production structure 
function for neutrino induced scattering has been computed in x space in Refs. [611162] : 



+ 2e- 



a. 



r f dz 



27r \ J z 



(9) 



where 



"'A , _ Q 



In Eq. Q, C^g^'^ includes the contributions in which the gluon splits into a s and a c 

quark, both of which contribute to at NLO. The Feynman diagrams for the LO 

and NLO gluon-induced subprocesses are shown in Figs. [H [H 

The x-space expressions for the O (as) charged current coefficient functions in Eq. ([9]) 
are given in Refs. [61..62J . The quark coefficient function can be separated into a delta 



^Partial knowledge of 0{ag) massive terms is available in the form of the asymptotic FgS] and threshold 
limits .60.. 
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c 



Figure 5: Feynman diagram for tire LO contributioir to F2^"^')'C;c FFm scheme. Thiclt solid 

liires indicate a heavy quark (charm) and thin solid lines a light quark (strange). 





Figure 6: Feynman diagrams for the NLO (O (as)) gluon-induced contribution to JT'^'"')'^^ in the 
FFNS scheme. 

function piece, a regular piece and a singular piece regulated with the usual plus prescrip- 
tion, 



Cif' {z, Q^ A) = (A) 6{l-z) + C^^ (A, z) + C^^^ (A, Q^) 



(11) 



The explicit expressions for the different pieces are the following. For the delta term we 
have 

1 1 + ^ 

2A + T + ^^^ 



Ka = (l-A)ln(l-A) /A . 
The regular piece can be written as 



(12) 
(13) 



cft)(A,z) = Cf 



(1 + z) (2 In {I- z)- In (1 - \z)) 



+ ( 2z + 2 - - 

z 



2 \ 1 

--1-Z 

z / 1 — Az 



[l + z^) Inz 
1 - z 



(14) 
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and finally the singular piece reads 

2 

+ 



1 + Q'^ + ml 21n(l-2) -ln(l- Az) 
1 I- z 



l-z 



i-z 2{i-\zy 



(15) 



where the first term is the contribution that depends on the factorization scale and is pro- 
portional to the qq splitting function. Separating the massive quark coefficient functions 
into the various contributions is important to properly evaluate their Mellin transforms, 
as will be discussed below. 

Finally, we give the expression for the FFN gluon coefficient function. In this case 
there are no singular terms and it reads 



a 



2,3 



{Z,Q\X) 



In — h In 



+ Tf[z^ + {l-zf^ (21n(l-z)-ln(l-Az)-lnz) 
+ f8-18(l-A) + 12(l + A)2)z(l-z)+ (^i^- 
1 - \z 



+ (l-A)zln 



(l-A)z 



(6A - I2\^z) 



(16) 



Again the last term in the first line is the scale-dependent contribution and is propor- 
tional to -Pgg'*. Note that both the diagrams shown in Fig. [6] contribute [61] . Analogous 
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expressions for the charged current F^^c. and Fl^c structure functions can be found in 
Refs. [611162]. 

As in the case of neutral currents, the massless limits of the FFN structure functions 
is easily obtained. For the massive -^2"^ structure function it has the structure 



F. 



(n,,0),CC 



2,c 



(x,Q2) 



2xs (x, Q 



+ 2x 



+ a 



as (Q2) 



27r 

(ni,0),l 
2,3 



dz 

z 



a 



{z,Q',X)g{-,Q' 



where 



a 



2Ji 



\z^Q\X) = Ctf'5{l-z) + c\:f\z) + 



(ni,0) 



(^) 



h,5 



c, 



(nifi) 



» = Cf 



(1 + z) In (l-z) - -^^ — + 3 + 2z 



l-z 



(17) 

(18) 
(19) 
(20) 
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h,s 



z) = Cf 



ln(l - z)\ 3 / 1 



1 - z 



2\l-z 



(21) 



and for the gluon 



2Tf 



[z' + (1 - zf) In 



1 - z 



+ 8z(l-2)-l 



+ Tf{z^ + {l-zf)\n^. 



(22) 



For completeness, we provide also the ZM-VFN quark coefficient functions for quarks 
and gluons, 



-(1 + z) In(l-z) 



ln(l - z] 
l-z 



(1 + z'^)\nz 



1 



2 \1- z 



+ 3 + 2z + (5 (1 - , 
l-z ^ ' \ 3 2 



vr^ 9 



a 



2,9 



{z) = Tf 



z^ + {1- zy) In 



+ (8z(l 



1) 



(23) 



(24) 



Note that the above gluon coefficient function, Eq. ()24p . is defined according to the nota- 
tion of Ref. [63], that is, it corresponds to the production of a single quark or antiquark. 

Comparing the FFNSO and ZM-VFN coefficient functions we find that for the gluon 
piece the following relation holds 



a 



{ni,0),l 
2,9 



2a 



2,9 ' 



z)+Tf(z^ + il-zf)ln9l, 



(25) 



where the overall factor 2 is due to the fact that the ZM coefficient function, Eq. ()24p . has 
been defined for a single quark, while in Eq. (jl6p the gluon coefficient function accounts 
for the production of two quarks {s and c). Note also the presence of the usual collinear 
logarithm. For the quark piece we find 



a 



2,h 



(z) = a 



2,h ' 



(26) 



without any collinear logarithm. 

The definition of the heavy CC structure function in the ZM scheme is not unique: here 
we define it as the contribution to the structure function which includes all contributions 
to the inclusive structure function which survive when all CKM elements but |V^s| are set 
to zero. With this definition, both the leading-order processes cW~^ — t- s and cW~^ — )■ s 
contribute to it (see Fig. [7]). This definition coincides with the experimental one because 
the struck charm antiquark must be accompanied by an (observed) charm quark, and it 
is free of mass singularities. The gluon initiated NLO contributions remain those shown 
in Fig. O The structure function in the massless scheme above charm threshold is then 
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Figure 7: Feynman diagrams that contribute to F. 



(ni+l),CC 



in the ZM-VFN scheme at leading 



order. The NLO diagrams are the same as in the FENS scheme. 



given by 

^^(n,+i),cC(^^^2) ^ 2x{s{x,Q^)+c{x,Q^)) 

a.^ (Q^) ( f dz 



+ 2x- 



2tt 



X 



z 



(27) 



The ZM-VFN massless coefficient functions have been defined in Eqs. ()23p -()24 p . Note the 
factor two in front of the gluon coefficient function, to account for the production of two 
quarks in the two NLO subprocesses of Fig. [6j 

Finally, the various schemes can be combined to construct the FONLL-A structure 
functions. As in the NC case, we define the FONLL structure function as follows 

Fi™)'--(x,Q^) = + e{Q^- ml) (l - ^) Q^) (28) 

f(J'^^(x, Q2) = f£'+^)'^^(x, Q2) - Fi;,"°)'^^(:r, Q') , (29) 

where as in the case of neutral currents we use the damping factor as default threshold 
prescription. 

Using the explicit expressions derived in the previous section for the difference between 
the ZM and FFNSO coefficient functions, Eqs. ()26p and ()25p . we can write the difference 
term as 

Fi;^,)'^^ = 2xc- (x, Q') - 2x^ In ^ ^T^ {z' + (1 - zf) g Q^) + O(a^) , (30) 

where we have used the fact that the heavy quark distribution is 0{as)- Now, it is easy 
to see explicitly that, in the region where L = In Q^/m^ is not large, the "difference" 
term is of order 0{a1): to first order in ag the FONLL expression coincides with the 
massive-scheme one also for charged current scattering. The use of the leading-order 
QCD evolution equations immediately leads to 

c(x, Q2) = c-(x, Q^) = In ^ ^r;(z2 + (1 - zf)g f-, Q^) + 0{al) . (31) 



27r mi .L z \z 
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Inserting this expansion in Eq. (j30p . it is trivial to check the exphcit cancellation of the 
0{as) terms, that is, that near the heavy quark threshold the difference term is of order 



(d),CC 



2.C 



The final FONLL-A expressions for the charged current charm production structure 
function is given by 



+ 2^ 



^(FONLL),CC(^^^: 
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(32) 
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1 _ !!k 



2x 



as {Q^) dz 



27r 



/ -Tf{z' + {l-zf)g(^,Q'^ 

Jx Z ' \Z 



It can be easily verified that Eq. (j32|l reduces to the FEN scheme, Eq. ([9]) at the heavy 
quark threshold = m^, and to the ZM-VEN expression Eq. (p7|) in the asymptotic 
region Q"^ » m^. 

The above derivation generalizes straightforwardly to the other relevant charged cur- 
rent structure functions xF^^ and F^^, as well as to the case with a general CKM quark 
mixing matrix. Note that in all the results shown below the standard CKM mixing has 
been assumed, with the CKM matrix elements set to their PDG values |64j . 

Now that we have defined the EONLL-A general mass scheme for charged current 
structure functions, we can compare the various schemes (ZM, EENS, EONLL-A) in the 
kinematic region that is most relevant in the global PDE analysis, namely the region 
covered by the NuTeV dimuon measurements |65j (see Eig.[T]). In Eig.[8]we show the results 
of such a comparison between various schemes for charm production in neutrino-induced 
charged current scattering. Results are compared at the level of the phenomenologically 
relevant charm production reduced cross section, defined as [9]: 



1 



dx dy 



{x,y,Q^ 



GIMn 



27r(l + QVM^ 



2 \2 



2M%x^y^ 



1 + ^ 



+ f]F!^^:\x,Q 



±y_xF3^f)(x,Q2) 



(34) 



with Q2 = 2MNE^xy and y± = 1 ± (1 - yf. In Eig. [8] we compare the various schemes in 
some representative bins of the NuTeV dimuon kinematics [65J. PDEs and other settings 
are those of the Les Houches heavy quark benchmark comparison [53j. We observe that in 
the kinematic region of neutrino data (both inclusive CHORUS data and dimuon NuTeV 
data), the EONLL-A result is very close to the EEN scheme computation, and it only 
begins to differ from it at the highest energies, where resummation of charm mass coUinear 
logarithms begins to become relevant. 



(33) 
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Figure 8: Comparison of different schemes for charm production in neutrino-induced DIS. The 
kinematic range is representative of the NuTeV dimuon data range. We compare the ZM-VFN, 
FFN and FONLL-A schemes at the level of the neutrino induced charm production cross section, 
Eq. The settings are the same as those of the Les Houches heavy quark benchmark comparison 

Even if the differences between the FFN and FONLL-A schemes for charged current 
scattering in the NuTeV kinematic region are moderate, as shown in Fig. [HI they become 
rather more important at small-x and medium-large Q^, where the charm and gluon PDFs 
become larger. To illustrate this, in Fig. [9] we compare the charged current charm structure 
function -^2^^ as a function of for two different values of x. Notice in particular that 
at very small x the FONLL-A expression is essentially the massless result. However, 
producing dimuons at x ~ 10^'^ and ^ 10 GeV^, where differences are larger, requires 
a fixed target neutrino experiment with a neutrino beam with energy in the multi-TeV 
range, which is not foreseen in the near future. Therefore one can conclude that any 
reasonable general-mass scheme for charged current scattering will be very close to the 
FFNS in the region of experimental data. 

The FONLL-A calculation of charged current structure functions has been imple- 
mented in a x-space code, FONLLdisCC, that we will use for benchmarking purposes. 
This is the analogue of the FONLLdis code for neutral currents [66], however is rather 
simpler since the unknown O (a^) massive coefficient functions do not have to be im- 
plemented. Our implementation of the FFNS calculations has been benchmarked with 
the corresponding results of the MSTW08 code [67], finding perfect agreement. We have 
also compared the FONLL-A and MSTW08 general-mass schemes for charged currents, 
finding qualitative agreement but some quantitative differences. A detailed comparison 
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LH Benchmark settings, x=0.033 



LH Benchmark settings, x=0.0001 




Figure 9: Comparison of different schemes for charm production in neutrino-induced DIS. We 
show the structure function in the massless, massive and FONLL-A schemes; in this case the 
FONLL- A expression is given by Eq. (15^ . The settings are the same as those of the Les Houches 
heavy quark benchmark comparison j53j . 

between different general-mass schemes for charged current structure functions, analogous 
to the Les Houches benchmarks for neutral current structure functions [53] is still missing 
and would be highly desirable. 
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2 

Xtot 


1.16 


{E) ± GE 


2.24 ± 0.09 


It? \ _i_ ^ 


z.z2 ± 0.11 




2.28 ±0.12 


(TL) ± C7TL 


(1.6 ±0.6) 10^ 




1.25 ±0.09 




11.3% 




4.4% 




0.18 




0.56 



Table 3: Table of statistical estimators for NNPDF2.1 with N^^^ = 1000 replicas. The total 
average uncertainty is given in percentage. All the and E values have been computed using the 
the same to covariance matrix [fiSJ used for minimization. 

4 Results 

In this Section we present the NNPDF2.1 parton determination. First, we discuss the 
statistical features of the fit, then we turn to a comparison of NNPDF2.1 PDFs and 
uncertainties with other PDF determinations and with previous NNPDF releases. A 
detailed comparison between the NNPDF2.1 and 2.0 sets follows, in which we discuss one 
by one the impact of the differences between the two fits, due both to the choice of dataset 
and to the different theoretical framework. The implications of the NNPDF2.1 set for 
LHC observables are discussed in the next sections. 

4.1 Statistical features 

Statistical estimators for the NNPDF2.1 fit are shown in Tab. [3] for the global fit and 
in Tab. H] for individual experiments. In Tab. H] the values for NNPDF2.0 are also 
shown for comparison. As in Ref. [lOJ, x^^^ is computed comparing the central (average) 
NNPDF2.1 fit to the original experimental data, (x'^^*^^) is computed comparing to the 
data each NNPDF2.1 replica and averaging over replicas, while {E) is the quantity which 
is minimized, i.e. it coincides with the computed comparing each NNPDF2.1 replica 
to the data replica it is fitted to, with the three values given corresponding to the total, 
training and validation data sets. It is important to observe that all values of shown 
in Tabs. [3ll3] are obtained using the covariance matrix with normalization uncertainties 
included according to the method of Ref. [68] (also given as Eq. (1) in Ref. [TO]). In 
Tabs. 9-10 of Ref. [TO] all values of Xtot (x^^'^^) were instead given with defined 
using the "standard" covariance matrix (given e.g. in Eq. (52) of that reference), which 
includes normalization uncertainties less accurately than the to covariance matrix. This 
was done in order to ease comparison between the results of Ref. [TO] and the NNPDFl.x 
PDFs, in which the to method was not yet used and normalization uncertainties were not 
fully accounted for. The values of Xtot NNPDF2.0 shown here in Tab. [3][3]has been 

recomputed using the to covariance matrix in order to ease comparison with NNPDF2.1 
and with other PDF determinations which also include normalization uncertainties albeit 
with various other methods. The value of Xtot fo'^ the NNPDF2.0 global fit computed using 
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Experiment 




2 

X2.0 


(E) 










NMC-pd 


0.97 


1.04 


2.04 


1.9% 


0.5% 


0.03 


0.37 


NMC 


1.73 


1.73 


2.79 


5.0% 


1.5% 


0.16 


0.71 


SLAG 


1.27 


1.42 


2.34 


4.4% 


1.6% 


0.31 


0.79 


BCDMS 


1.28 


1.30 


2.33 


5.7% 


2.3% 


0.47 


0.60 


HERAI-AV 


1.07 


1.15 


2.15 


2.5% 


1.2% 


0.06 


0.35 


CHORUS 


1.15 


1.24 


2.23 


15.1% 


4.7% 


0.08 


0.32 


FLH108 


1.37 


1.50 


2.36 


72.0% 


4.0% 


0.64 


0.67 


NTVDMN 


0.76 


0.73 


1.77 


21.1% 


14.1% 


0.04 


0.62 


ZEUS-H2 


1.29 


1.33 


2.32 


13.4% 


1.2% 


0.27 


0.51 


ZEUSF2C 


0.78 




1.80 


23.3% 


3.1% 


0.08 


0.41 


H1F2C 


1.50 




2.52 


17.3% 


3.0% 


0.30 


0.40 


DYE605 


0.84 


0.87 


1.92 


22.3% 


7.9% 


0.47 


0.76 


DYE866 


1.27 


1.29 


2.37 


20.1% 


9.2% 


0.20 


0.52 


CDFWASY 


1.86 


1.84 


3.08 


6.0% 


4.4% 


0.51 


0.75 


CDFZRAP 


1.65 


1.85 


2.80 


11.5% 


3.6% 


0.82 


0.72 


DOZRAP 


0.60 


0.60 


1.62 


10.2% 


3.1% 


0.53 


0.76 


CDFR2KT 


0.97 


1.01 


2.10 


22.2% 


4.0% 


0.78 


0.57 


D0R2CON 


0.84 


0.86 


1.92 


16.8% 


4.5% 


0.77 


0.59 



Table 4: Same as Table |3] for individual experiments. All estimators have been obtained with 
A^rep = 1000 replicas. Note that experimental uncertainties are always given in percentage. In the 
second and third column the NNPDF2.1 and NNPDF2.0 set [10] have been computed with the 
to prescripion. 



the to method, to be compared to the NNPDF2.1 value of Tab. [3l is Xtot = 1-23 (very 
close to the the value Xtot — l-^l of Tab. 9 in Ref. [lOj, computed with the "standard" 
covariance matrix). 

The NNPDF2.1 PDF fit has the following noticeable features: 

• The quality of the global fit as measured by the value = 1-16 is rather better 
than for the NNPDF2.0 fit without heavy quark mass effects. 

• As compared to the NNPDF2.0 results, the quality of the fit to all datasets im- 
proves or remains similar. The most noticeable improvements can be found for the 
HERA-I average dataset and for CHORUS. The improvement in the description of 
HERA data arises both from the improved heavy fiavour treatment and the more 
conservative kinematical cuts. 

• An excellent description of the combined HERA-I inclusive data, = 1.07, is 
obtained. Similarly, a reasonable description of the HERA charm structure function 
data is achieved. 

• The quality of the fit to hadronic data is not affected by the use of the FONLL-A 
GM scheme for deep-inelastic observables, as it can be seen by comparing the second 
and third column of Table [H 

The distribution of x^^'^K E^'' and training lengths among the A^'rep = 1000 NNPDF2.1 
replicas are shown in Fig. [TU] and Fig. [TT] respectively. While most of the replicas fulfill 
the stopping criterion, a fraction (~ 12%) of them stops at the maximum training length 
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y^''*' distribution 



E,, distribution 




Figure 10: Distribution of ^^^(fe) (left) and E^^' (right), over the sample of N^-ep — 1000 repUcas. 



which has been introduced in order to avoid unacceptably long fits. This causes 
some loss of accuracy for outliers fits (i.e. those in the tail of the distribution): we have 
checked that as is raised more and more of these replicas stop, and that the loss of 

accuracy due to this choice of value of N^^^ is reasonably small, in that the features of 
the global fit change very little if is raised. 

^Distribution of training lenghts | 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 

^00 10000 15000 20000 25000 30000 

Training lenght [GA generations] 

Figure 11: Distribution of training lengths over the sample of N^-ep = 1000 replicas. 




It is instructive to compare the quality of the fit with the corresponding results obtained 
in the recent CTIO analysisjfl In Table [5] we compare the of the common sets in 
NNPDF2.1 and CTIO, along with the number of data points in each fit (which differ 
because of different kinematic cuts, see Tabled]). It should be borne in mind that the 

is defined in a somewhat different way by the CTEQ/CT group, specifically, but not 
only, in what concerns the treatment of normalization errors (see Ref. [16j): hence this 
comparison should be taken with care. From this comparison, we can see that the two 
sets have a comparable fit quality to fixed target DIS, CTIO being somewhat better for 

®We thank Pavel Nadolsky for providing us with these numbers. 
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NNPDF2.1 


CTIO 


Experiment 


A^dat 




A^dat 




NMC-pd 


132 


0.97 


121 


1.28 


NMC 


221 


1.73 


196 


1.71 


BCDMSp 


333 


1.28 


337 


1.14 


BCDMSd 


248 


1.15 


250 


1.12 


HERAI-AV 


592 


1.07 


579 


1.17 


NTVnuDMN 


41 


0.50 


38 


0.94 


NTVnbDMN 


38 


0.42 


33 


0.91 


DYE605 


119 


0.85 


119 


0.81 


DYE866p 


184 


1.31 


184 


1.21 


DYE866r 


15 


0.77 


15 


0.64 


CDFZRAT 


29 


1.62 


29 


1.44 


DOZRAP 


28 


0.59 


28 


0.54 


CDFR2KT 


76 


0.97 


76 


1.55 


D0R2CON 


110 


0.84 


110 


1.13 



Table 5: Comparison of per data point for experiments which are common to the NNPDF2.1 
and CTIO PDF determinations. For each PDF set the number of data points obtained with the 
kinematic cuts of Table [T] is given. 

BCDMS proton and NNPDF2.1 rather better for NMC deuteron/proton ratio. The fit 
to HERA-I and Tevatron jet data is rather better in NNPDF2.1. Comparable fit quality 
to the Drell-Yan and vector boson production data is obtained in the two cases, with 
somewhat smaller in the CTIO fit. No comparison is attempted for the HERA F| data 
because of the very different kinematic cuts used in the two fits. A similar comparison to 
MSTW08 would be less significant because in the MSTW08 fit correlated systematics are 
not included in the covariance matrix for some datasets. 

4.2 Parton distributions 

The NNPDF2.1 PDFs are compared to the previous NNPDF2.0 PDFs in Figs. [I2] (singlet 
sector) and [13] (non-singlet sector) . 

• The singlet PDF at medium and small-x is rather similar in the two cases, but it is 
somewhat larger in the NNPDF2.1 set. 

• Thanks to the new positivity constraint on F2 in NNPDF2.1, the gluon remains 
always positive even at the largest values of x, where occasionally went very slightly 
negative in NNPDF2.0. 

• In NNPDF2.1 the small-x gluon is larger than in NNPDF2.0. We wih show that this 
arises from the use of a GM scheme as compared to the ZM scheme in NNPDF2.0. 
Also, the medium and small-x gluon has a somewhat larger uncertainty in NNPDF2.1 
as compared to NNPDF2.0. We will show below that this uncertainty increase is 
due to the new kinematic cut. 
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Figure 12: Comparison of NNPDF2.1 and NNPDF2.0 singlet sector PDFs, computed using iVrop 
1000 replicas from both sets. All error bands shown correspond to one-cr. 



As expected differences in the large-x valence PDFs are rather modest. The GM 
scheme does not affect the large-x PDFs but the cross-talk induced by the sum rules 
and other constraints induces small modifications also in the valence sector, always 
well below the one-o" level. 

The strange PDF in NNPDF2.1 is somewhat smaller than in NNPDF2.0. This may 
appear surprizing as the main expected effect of the inclusion of heavy quark mass 
effects is a suppression of charm which then leads to an enhancement of all other 
PDFs. However, in NNPDF2.0 for dimuon data (which have a sizable impact on 
strangeness) instead of a pure zero-mass scheme, the so-called improved zero-mass 
(IZM) approximation of Ref. [SI] was used to approximate heavy quark mass effects. 
It turns out that this IZM method actually overestimates heavy quark mass effects, 
thus leading to a slight over-suppression of strangeness in NNPDF2.0. We will check 
explicitly below (see Fig. [8]) that when comparing NNPDF2.1 to a pure zero-mass 
fit strangeness is somewhat enhanced as one would expect. 

We see from Fig. [13] that the strange asymmetry in NNPDF2.1 is very close to that 
of NNPDF2.0. An important result of the NNPDF2.0 analysis was that the strange 
asymmetry s~{x, Q^) was of the proper size to completely cancel the so-called NuTeV 
anomaly with rather reduced uncertainties. It is clear that this holds true also with 
the updated NNPDF2.1 set, confirming the results of the analysis of [9j that showed 



28 




Figure 13: Same as Fig. [T2] for the non-singlet sector PDFs. 



that heavy quark mass effects have a very moderate impact in the determination 
of the strangeness asymmetry. The imphcations for the NuTeV anomaly will be 
discussed in Sect. 15.41 

The comparison between NNPDF2.0 and NNPDF2.1 is further quantified by the com- 
putation of the distance PH] between the two sets, shown in Fig. UM Note that d ~ 1 
corresponds to two sets of replicas which come from the same underlying probability dis- 
tribution, while (using A^rep = 100 replicas) d ~ 7 corresponds to a one-o" difference (see 
Appendix A of Ref. [lOj). One concludes that while clearly the two sets do not come from 
the same underlying distributions, all PDFs but the strange are consistent at the one-cr, 
and even the strangeness is consistent at the 90% confidence level. The largest differences 
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Figure 14: Distance between the NNPDF2.0 and NNPDF2.1 parton sets. All distances are 
computed from sets of N^cp = 100 replicas. 



are seen in the medium-x strangeness and to a lesser extent in the medium and small-x 
gluon. 

The differences between NNPDF2.0 and NNPDF2.1 PDFs at the initial scale displayed 
in Figs. [T2lfHl once evolved up to the W and Z scale, are sufficient to lead to differences 
between gluon and light sea quark distributions up to the one-cr level at small-x, as shown 
in Fig. [T5]where we plot the NNPDF2.0/NNPDF2.1 ratio for individual light flavours and 
the gluon at = 10^ GeV^. 

The NNPDF2.1 PDFs are compared to the other global PDF sets CTIO [l6] and 
MSTW08 [21] in Figs. [T6]|17l which is interesting to contrast to the analogous plot which in 
Ref. [lOj (Figs. 18-19 of that reference) compared the NNPDF2.0, CTEQ6.6 and MSTW08 
PDF sets. 

• The general agreement of the gluon in the medium-/small-x region is improved, both 
because the central value of NNPDF2.1 is now in better agreement with MSTW08 
(most likely due to the inclusion of heavy quark mass effects) and also because the 
CTIO central value and especially uncertainty are in much better agreement with 
the wider NNPDF and MSTW uncertainties, due to the use of a more flexible gluon 
parametrization in CTIO with respect to CTEQ6.6. The large-x gluon however is in 
marginal agreement. 

• The small changes in valence and triplet distributions between NNPDF2.0 and 
NNPDF2.1 go anyway in the direction of improving the agreement with the other 
global sets. 
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Figure 15: Comparison between NNPDF2.1 and NNPDF2.0 liglit quark and gluon PDFs at 
= lO'* GeV^. All curves are shown as ratios to the central NNPDF2.1 result. 



• The strange PDFs are quite different, presumably due to the fact that a much less 
flexible parametrization is adopted by CT/CTEQ and MSTW in comparison to 
NNPDF. 

• The medium-x singlet is in marginal agreement. This may be a byproduct of the 
poor agreement in strangeness. 

The effect on LHC observables will be discussed in Sect. \5\ where we will show that even 
though there is generally a reasonable agreement between global sets, there remain some 
significant differences, mostly related to the rather different large-x gluon in CTIO as 
shown in Fig. [TBI 

4.3 Detailed comparison to NNPDF2.0: theoretical framework and dataset 

We now assess the separate impact of the inclusion of heavy quark mass effects and of the 
charm structure function data on the NNPDF2.1 PDF determination. 

First, we must discuss the impact of raising the kinematic cut Q^j^ within the ZM 
scheme; then we compare the ZM and GM fits; next we investigate the impact of including 
the HERA charm structure function data into the NNPDF2.1 analysis; finally we estimate 
the impact of ambiguities related to the treatment of heavy quarks. In each case, we will 
show the distances between PDFs as well as the PDFs that are most affected by each step. 
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• New kinematical cut Q^i^ 

In Fig. [18] we show the distance between NNPDF2.0 PDFs and a fit with the same 
dataset but with the new cut Q^^^ = 3 GeV^, denoted by NNPDF2.0RED (reduced). 
Also, in order to ease the subsequent discussion of the impact of heavy quark mass 
effects, in NNPDF2.0RED a pure ZM scheme is used for all observables, rather than 
the IZM scheme [5l] used for dimuon data in Ref. |10j . 

The largest distances correspond to the medium- strange PDFs and the small-x 
gluon. These PDFs are shown in Fig. [TUl The strange is rather smaller in the ZM 
as compared to the IZM scheme, where it was enhanced due to the approximate 
inclusion of charm suppression. The gluon is somewhat smaller at small-x and with 
rather larger uncertainties, due to the reduction in dataset at small- a; caused by the 
new kinematic cut. From the distances we see that the singlet is also modified, but 
one can check that this is completely due to the strange contribution to it. 

• Impact of the general-mass scheme 

The impact of the FONLL-A GM is assessed by now comparing the NNPDF2.0RED 
fit to a NNPDF2.1 fit without F2 data: the dataset is identical, the only difference 
is in the treatment of heavy quark masses. The distances between these two sets are 
shown in Fig. [20l 




Figure 16: The NNPDF2.1 singlet sector PDFs, compared with the CTIO and MSTW08 PDFs. 
The results for NNPDF2.1 have been obtained with A^rcp = 1000 replicas. All PDF errors are given 
as one-fj uncertainties. 



32 




■ NNPDF2.1 
I ICTIO 
r~1 MSTW08 



o 
— 0.1 





■ NNPDF2.1 
1 ICTIO 
f 1 MSTW08 






, , , 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 , . , , 1 , . , , 1 , . , , 1 , , , , 1 , , , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



■ NNPDF2.1 
I ICTIO 
□ MSTW08 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 







Wk NNPDF2.1 






1 ICTIO 






W\ MSTW08 













^NNPDF2.1 

I ICTIO 
C3 MSTW08 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 17: Same as Fig. [Tnifor tlie non-singlet sector PDFs. 



The impact of the inclusion of heavy quark masses is mostly on the small-x gluon 
and, to a lesser extent, on medium-x strangeness. These two PDFs are shown in 
Fig. [2TJ The GM scheme leads to a larger gluon for x < 2 • 10^^, as well as to a 
somewhat larger strangeness, but it leaves the singlet unaffected. This shows that 
indeed, as argued in Sect. 14.21 the relatively large total strangeness in NNPDF2.0 
was due to the use of the IZM approximation for dimuon data, which overestimates 
charm mass effects. 

• Impact of HERA F| data 

The impact of HERA data is estimated comparing the results of the NNPDF2.1 
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Figure 18: Distance between the NNPDF2.0 PDF set and a fit to the same data but with Qcut = 3 
GeV^ and the ZM-VFN scheme for all observables (NNPDF2.0 RED). All distances are computed 
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Figure 19: Comparison of the small-a; total strangeness and gluon in NNPDF2.0 and in 
NNPDF2.0RED (the distances are shown in Fig.lT^. 



fit with and without these data. The distances displayed in Fig. 1221 show that current 
F2 data have little effect: the two fits are almost statistically equivalent, with most 
distances of order one. This is partly due to the relatively large uncertainties on 
current these F2, and also to the fact that low x and data, which are most 
sensitive to the gluon PDF, are excluded by our kinematic cuts. Inclusion of 0{a'^) 
heavy quark mass effects (e.g. by means of the FONLL-B scheme) is necessary in 
order to take advantage of these data. 
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Figure 20: Distance between the NNPDF2.1 PDF sets in the GM and in the ZM schemes, in both 
cases without HERA data. All distances are computed from sets of iVj-cp = 100 replicas. 




Figure 21: Comparison of the small-a; total strangeness and gluon in NNPDF2.0RED and 
NNPDF2.1 without data (distances are shown in Fig. BUI . 



• Impact of threshold prescription in the GM scheme 

Finally, we have repeated the NNPDF2.1 fit with a different treatment of subleading 
terms in the inclusion of heavy quark mass effects: namely, we have used FONLL-A 
but without the threshold damping factor in Eqs. (|l|32p . Indeed, the benchmarking 
exercise of Ref. |53j suggest that the difference between these cases should provide 
a reasonable estimate of the spread of results obtained by including heavy quark 
masses according to different prescriptions. Distances are shown in Fig. [23J the 
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Figure 22: Distance between NNPDF2.1 PDF sets with and without HERA F^{x,Q'^) data. 
Distances have been computed from sets of iVrop = 100 replicas. 



PDFs that are most affected are the singlet and gluon PDFs at medium-x, shown in 
Fig. [231 Without damping factor, the F2 structure function is closer to the massless 
result even at moderate Q^, and this explains why the singlet PDF is somewhat 
smaller at medium-x. 
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Figure 23: Distance between the NNPDF2.1 reference set and tlie same set obtained witliout 
threshold damping factor in the computation of the FONLL-A structure functions. . 
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a{W+)Bi, [nb] 


a{W-)Bi^ [nb] 


a{Z")Bu [nb] 


NNPDF2.0 


5.84 ±0.14 


3.97 ±0.09 


0.91 ±0.02 


NNPDF2.1 


5.99 ±0.14 


4.09 ± 0.09 


0.93 ±0.02 


CTIO - as = 0.118 
CTIO - as = 0.119 


6.00 ±0.13 
6.04 ±0.13 


4.10 ±0.09 
4.13 ±0.09 


0.94 ± 0.02 
0.95 ±0.02 


MSTW08 - as = 0.119 
MSTW08 - as = 0.120 


5.91 ±0.11 
5.95 ±0.11 


4.16 ±0.08 
4.19 ±0.08 


0.94 ± 0.02 
0.95 ±0.02 





a{tt) [pb] 


a{H,m„ = 120 GeV) [pb] 


NNPDF2.0 


168 ±7 


11.59 ±0.22 


NNPDF2.1 


170 ±5 


11.64 ±0.17 


CTIO - as = 0.118 
CTIO - as = 0.119 


158 ±7 
161 ±7 


10.99 ±0.21 
11.17 ±0.21 


MSTW08 - as = 0.119 
MSTW08 - as = 0.120 


164 ±5 

168 ±5 


11.48 ±0.18 
11.69 ±0.18 



Table 6: Cross-sections for W, Z, tt and Higgs production at the LHC at y/s — 7 TeV and the 
associated PDF uncertainties. All quantities have been computed at NLO using MCFM for the 
NNPDF2.1, NNPDF2.0, CTIO and MSTW08 PDF sets. All uncertainties shown are one-a. 

5 Phenomenological implications 

In this Section we discuss the implications of the NNPDF2.1 set for LHC physics. We begin 
comparing the prediction for LHC benchmark cross sections obtained using NNPDF2.1 
and NNPDF2.0. This comparison allows us to assess the impact of heavy quark mass 
effects. We also compare to predictions obtained using the CTIO and MSTW08 sets, both 
using their preferred values of as and with a common value. We then compare parton lu- 
minosities relevant for LHC processes and determine correlations between PDFs and some 
observables. Next, we determine the correlation between PDFs and and discuss PDF 
sets with varying as, which are needed to compute the combined PDF+a^ uncertainties, 
and present sets in which PDF and Og uncertainties are pre-combined. Finally, we briefly 
revisit implications of the strangeness asymmetry on the NuTeV anomaly, confirming our 
previous result that the anomaly disappears once the strangeness asymmetry is properly 
determined. We conclude with a comparison of NNPDF2.1 results with published HERA 
Fl and F2 data and predictions for the upcoming and F2 HERA combined datasets. 

5.1 LHC benchmark cross-sections 

The assessment of the theoretical uncertainties on LHC standard candles is especially 
important now that the first 7 TeV LHC results on inclusive cross-sections are appear- 
ing In this Section we present results at ^/s =7 TeV and ^/s =14 TeV for W^, 
Z^, ti and Higgs production via gluon fusion with mu = 120 GeV. All observables are 
computed at NLO QCD using MCFM ^UM- 

In Tables [6] and [71 and the corresponding Figs. [25] and [271 we compare the predictions 
for these cross-sections obtained using the NNPDF2.1, NNPDF2.0, CTIO and MSTW08 
sets. In the case of the last two sets, we show results both using the respective default 
value of as{Mz) and at the common value of as{Mz) = 0.119, obtained using the PDF 
sets of Refs. [75l[76]. 

Predictions obtained using the NNPDF2.0 and NNPDF2.1 PDF sets mostly differ 
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Figure 25: Graphical representation of the results of Table IHl 



because of heavy quark mass effects, but other differences such as different kinematic 
cuts also play a role. As can be seen from Fig. [25] and Table El the differences between 
NNPDF2.0 and NNPDF2.1 are at most at the one-cr level for and Z production 7 
TeV, while predictions for the ti and Higgs are essentially unchanged: these observables 
are only minimally affected by the heavy quark treatment. 

NNPDF2.1 predictions are in rather good agreement with MSTW08 for all observables, 
though differences with CTIO are somewhat larger, especially for observables which are 
most sensitive to the gluon distribution, like Higgs and production. The use of a common 
value for the strong coupling ag leads to better agreement between predictions, especially 
for processes which depend on ag already at leading order such as Higgs production in 
gluon fusion [18j. In Fig. [26] first measurements by the ATLAS and CMS experiments [691 
[70] are also shown: with their large uncertainties, dominated by the current large (0(11%)) 
luminosity uncertainty, they cannot yet provide constraint on PDFs. 

At -y/s =14 TeV we expect the effect of the heavy quark treatment to be larger: results 
are collected in Fig. [27] and Table [7] In this case, the upwards shift in the the and Z 
cross-sections from NNPDF2.0 to NNPDF2.1 is at or just above the one-a level, while as 
before Higgs and top-pair cross-sections are essentially unchanged. The comparison with 
CTIO and MSTW08 is similar as before, but with the agreement somewhat better for the 
Higgs and somewhat worse for top. 

Related important observables at the LHC are the W~^/W~ and W/Z cross-section 
ratios. These have generally reduced experimental uncertainties, since e.g. normalization 
uncertainties cancel in the ratio. Predictions at 7 and 14 TeV for NNPDF2.1, CTIO 
and MSTW08 are compared in Fig. [28] For these observables the dependence on as is 
negligible. The agreement for cross-section ratios seems to be worse than for total cross- 
sections: for example for the W'^ /W~ ratio CTIO and NNPDF2.1 are in good agreement 
but MSTW08 is lower by more than two-cr. The agreement for the W/Z ratio is better 
but still marginal at 7 TeV. 
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Figure 26: LHC measurements of the W^,Z and tt cross-sections at ^/s —7 TeV from the 
ATLAS [51] and CMS experiments [TD] compared to the predictions of Fig. [55] 



The correlation between PDFs and physical observables quantifies the relevance of 
each PDF for different observables (and conversely) as a function of x [H[I5]- As an 
illustration with NNPDF2.1, we have computed the correlations between PDFs and vector 
boson production total cross-sections and their ratios at LHC 7 TeV. Results are shown 
in Fig. [29] The total W and Z cross-sections are as expected mostly correlated with the u 
and d sea quarks and anticorrelated with the strange quarks, the correlation with the gluon 
(and the heavy flavours generated dynamically from it) being milder. It is also interesting 
to note that correlations between PDFs and the corresponding physical observable are 
only moderately reduced in the W/Z ratio as compared to the individual cross-sections, 
and they are instead almost suppressed in the /W~ ratio. This observation suggest 
that the latter ratio should be less sensitive to PDF uncertainties. 

As for NNPDF2.0 [lO], we have produced variants of the NNPDF2.1 fit based on 
reduced datasets: DIS only, DIS and inclusive jet data only. These fits are useful to study 
the impact on PDF of the various observables used in the global fit. Results for LHC 
cross-sections at 7 TeV determined using these fits are collected in Table [5] For these very 
inclusive observables, it turns out that a purely DIS fit already provides a rather good 
approximation, though this need not be always the case for other observables. 

The inclusion of heavy quark mass effects has a theoretical ambiguity due to subleading 
terms. A full study of theoretical uncertainties on PDFs has never been performed and 
goes beyond the scope of this paper. However, we provide here a first estimate of the 
uncertainty related to the inclusion of heavy quark mass effects to 0{as) by comparing 
results obtained from the three sets discussed in Sect. 14.31 NNPDF2.0 RED (without heavy 
quark mass terms, but the same kinematic cuts NNPDF2.1), the default NNPDF2.1, and 
NNPDF2.1 without damping terms in the FONLL-A method. Results are shown in Table[2] 



40 



LHCUTeV 



LHC UTeV 



LHC UTeV 




Figure 27: Graphical representation of the results of Table [71 



and in Fig. [301 As expected, results obtained without damping prescription sit half way 
between NNPDF2.0RED and NNPDF2.1 at 7 TeV, and closer to the latter at 14 TeV. 
As discussed in Sect. 14.31 the difference between the NNPDF2.1 results with and without 
damping terms can be taken as a conservative estimate of the theoretical uncertainty 
associated to the uncertainty in the inclusion of heavy quark mass effects to 0{as)- 

5.2 Parton luminosities 

The processes discussed in Sect. 15.11 are a small subset of the LHC observables which are 
sensitive to PDFs. A detailed, systematic study of these would be quite interesting; how- 
ever, a good deal of information can be gathered by simply studying parton luminosities. 
Following Ref. [77], we define the parton luminosity 

ci>.,- (Mi) = - [' (xi, A4) /, (r/:ri. Mi) , (35) 

where fi{x,M'^) is a PDF and r = M'j^/s. We consider in particular the gluon-gluon 
luminosity, the various heavy quark-antiquark luminosity, and the quark-gluon and quark- 
quark luminosity respectively defined as 

^,9^Y.^m-, ^m^Y.^m.- (36) 

i=l i=l 

Parton luminosities form NNPDF2.1, NNPDF2.0, CTIO and MSTW08 sets, normal- 
ized to the NNPDF2.1 central prediction, are shown in Fig. [31] for LHC y/s = 7 TeV, all 
determined with ag (Mz) = 0.119. In Fig. [32] we also compare directly the relative uncer- 
tainties on the luminosity for each set. These comparisons show good agreement between 
global fits at the one-u level, although in some cases, such as the gluon-gluon luminosity 
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a{W+)Bi, [nb] 


a{W-)Bi^ [nb] 


a{Z")Bu [nb] 


NNPDF2.0 
NNPDF2.1 


11.59 ±0.27 
12.00 ±0.27 


8.56 ±0.17 
8.84 ±0.17 


1.94 ±0.04 
1.99 ±0.04 


CTIO - Q, = 0.118 
CTIO - = 0.119 


12.20 ±0.30 
12.31 ±0.30 


9.00 ±0.22 
9.07 ±0.22 


2.03 ± 0.05 
2.05 ± 0.05 


MSTW08 - a, = 0.119 
MSTW08 - as = 0.120 


11.95 ±0.22 
12.06 ±0.22 


9.03 ±0.17 
9.10 ±0.17 


2.01 ±0.04 
2.03 ±0.04 





a{ti) [pb] 


a{H,mH = 120 GeV) [pb] 


NNPDF2.0 
NNPDF2.1 


942 ± 21 
946 ± 19 


37.3 ±0.50 
37.5 ±0.40 


CTIO - as = 0.118 
CTIO - as = 0.119 


880 ± 21 
895 ± 21 


36.32 ±0.80 
36.90 ± 0.80 


MSTW08 - as = 0.119 
MSTW08 - as = 0.120 


917 ± 18 
934 ± 18 


37.78 ± 0.50 
38.43 ± 0.50 



Table 7: Same as Table E] for the LHC at ^ 14 TeV. 





a{W+)Bi,, [nb] 


a{W-)Bi^ [nb] 


a{Z")Ba [nb] 


NNPDF2.1 DIS 
NNPDF2.1 DIS+JET 
NNPDF2.1 


6.03 ±0.11 
6.03 ±0.12 
5.99 ±0.14 


4.15 ±0.08 
4.14 ±0.08 
4.09 ± 0.09 


0.940 ±0.014 
0.939 ±0.015 
0.932 ± 0.020 





a{tt) [pb] 


a{H,m„ = 120 GeV) [pb] 


NNPDF2.1 DIS 
NNPDF2.1 DIS±JET 
NNPDF2.1 


167 ±7 
170 ±5 
170 ±5 


11.66 ± 0.21 
11.66 ± 0.22 
11.64 ± 0.17 



Table 8: Cross-sections for W, Z, tt and Higgs production at the LHC at ^/s = 7 TeV and the 
associated PDF uncertainties for the reference NNPDF2.1 set compared to those obtained using 
sets determined from reduced datasets: DIS only, DIS+JET. 

at intermediate invariant masses, the agreement is only marginal. Uncertainties blow up 
both at very small and large values of Mx for all sets. Differences between different sets 
are larger in these regions: for example NNPDF2.1 uncertainties at small Mx are rather 
larger, for luminosities that involve the gluon PDF. The peculiar behaviour of the bot- 
tom luminosity for NNPDF2.0 is due to the fact that in this set rrib = 4.3 GeV, while 
rrih = 4.75 GeV for all other sets. The dependence of results on the values of the heavy 
quark masses will be discussed in Sect. El below. 

5.3 The value of (M^) 

We now consider the correlation between NNPDF2.1 partons and the value of Og- To this 
purpose, we provide sets with ag (Mz) in the range from 0.114 to 0.124 in steps of 0.001. 
PDFs from fits performed using different values of ag are show in Fig. [33j Results are 
similar to those obtained with NNPDF2.0 [T7]: as expected, the most sensitive PDF is 
the gluon. To quantify this it is useful to compute the correlation between PDFs and 
(as defined in Eq. (82) of Ref. [IZ]). We determine it assuming the uncertainty on to 
be 5as = 0.0012 at the 68% C.L. Results are plotted in Fig. [M] as a function of x, both at 
= 2 GeV^ and = 10^ GeV^. Clearly, because of asymptotic freedom, correlations 
are weaker at high scale. 
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Figure 28: Comparison between predictions from different PDF sets for the /W and W/Z 
ratios at the LHC y/s =7 TeV (upper plots) and ^/s —14 TeV (lower plots). 



Following the procedure outlined in Ref . jl8j , it is possible to combine sets with different 
values of Og to compute the combined POF+a^ uncertainty on any given observable. This 
procedure has the advantage that both the central value and the uncertainty on the strong 
coupling are not fixed a priori but can be chosen by the PDF user. In order to simplify this 
procedure, we provide prepacked PDF sets with combined PDF+a^ uncertainty. Once a 
central value and uncertainty for as{Mz) are assumed, this is done by constructing a set 
of Nj-cp replicas, extracted from the original sets with different Og, in such a way that the 
prepacked sets contains a number of replicas for each value of Os which corresponds to a 
gaussian distribution with given mean and standard deviation |18J (of course, any other 
distribution could be used) . The statistical accuracy of the prediction obtained using the 
prepacked sets scales with the number of replicas A'rcp- 

We have produced prepacked PDF+a^ uncertainty sets with a<j (M^) = 0.119 and 
uncertainties Sa^ = 0.0012 and 5a^ = 0.002 as one-a errors. These values have been 
chosen to agree with the PDF4LHC recommendation [78] for the combination of PDF+as 
uncertainties. Sets with any other values are easily produced and are available upon 
request. For completeness, we have produced the same prepacked sets also for NNPDF2.0. 
We have checked that results for Higgs production in gluon fusion at LHC 7 TeV kinematics 
(which has rather large as uncertainties and correlations) become essentially independent 
of the number of replicas in the prepacked set provided A'rep > 100. For smaller number 
of replicas there is a certain loss of accuracy, so a minimum of A'^rep =100 is recommended. 

An important caveat in the usage of prepacked PDF sets is that some widely used 
codes, such as MCFM, assume the value of ag is the same for all PDFs in a given set. 
Prepacked sets cannot be used with these codes. Similar prepacked PDF sets could be 
prepared to include the uncertainty on other physical parameters, such as heavy quark 
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NNPDF2.1 , LHC7, Corr PDFs and o(W*) at Q"^ = M. 



NNPDF2.1 , LHC7, Corr PDFs and <j(Z) at = M^, 




1e-05 0.0001 0.001 0.01 



1e-05 0.0001 0.001 0.01 



NNPDF2.1 , LHC7, Corr PDFs and o(W*)/<j{W') at Q = M, 

1.5 I . . ^ 



NNPDF2.1 , LHC7, Corr PDFs and a(W)/a{Z) at Q = M„ 

1.5 I , , ^ 




1e-05 0.0001 0.001 0.01 



Figure 29: The correlation between PDFs and vector boson production total cross-sections (upper 
plots) and their ratios (lower plots) for LHC 7 TeV. Correlations for (not shown) are very 
similar. 



masses. An important limitation however is imposed by the current LHAPDF standard 
which assumes that all physical parameters except the the strong coupling take the same 
value for all PDFs in a given set. For this reason, only PDF+ag prepacked sets are 
provided for the time being. 

5.4 The NuTeV anomaly 

In previous NNPDF releases PlllUj we studied the implications that the determination of 
the strangeness asymmetry s~{x,Q'^) has on the so-called NuTeV anomaly [79]. These 
results are updated here. For the first moment of the strangeness asymmetry with 
NNPDF2.1 we find 

^ 2^ Iodxxs-{x,Q') ^ ^ [S-] ^ ^ 

f^dxx{u-{x,Q^) + d-{x,Q^)) [U-+D-] 

In Fig. [35] we show the NuTeV determination of the Weinberg angle [65] , uncorrected 
and then corrected for the strangeness asymmetry using the values from previous [milUj 
and the current NNPDF sets. The three corrected values are in excellent agreement with 
the electroweak fit and with each other, with the NNPDF2.0 and NNPDF2.1 values very 
close to each other, thereby showing that the impact of heavy quark mass effects on the 
determination of the strangeness asymmetry is very small. 
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7 TeV 


a{W+)Bi^ [nb] 


a{W-)Bi^ [nb] 


a{Z")Bii [nb] 


NNPDF2.1 


5.99 ±0.14 


4.09 ±0.09 


0.93 ±0.02 


NNPDF2.0 RED 
NNPDF2.1 FONLL-A plain 


5.81 ±0.13 
5.90 ±0.12 


3.98 ± 0.08 
4.03 ± 0.08 


0.91 ±0.02 
0.92 ±0.02 




14 TeV 


[nb] 


a{W-)Bi^ [nb] 


a{Z'')Bii [nb] 


NNPDF2.1 


12.00 ± 0.27 


8.84 ±0.17 


1.99 ±0.04 


NNPDF2.0 RED 
NNPDF2.1 FONLL-A plain 


11.57 ±0.25 
11.82 ±0.22 


8.57 ±0.17 
8.72 ±0.15 


1.93 ±0.04 
1.96 ±0.03 



Table 9: Cross-sections for W, Z, tt and Higgs production at the LHC at y/s = 7 TeV and the 
associated PDF uncertainties for the reference NNPDF2.1 set compared to those obtained using 
sets with different treatment of heavy quarks: NNPDF2.0RED, without heavy quark mass effects, 
and NNPDF2.1 FONLL-A plain with heavy quark mass effects but without threshold damping 
terms. 



LHC 7 TeV - Different HO si 



ZM-VFN FONLL-A-plam FONLL-A-damp 



LHC14TeV - Differenf HQ si 
ZM-VFN FONLL-A-plain FONLL-A-damp 



LHC 7 TeV - Diffetenf HQ sctiemes 



FQNLL-A-plain FQNLL-A-damp 



LHC 14 TeV - Different HQ scliemes 
ZM-VFN FONLL-A-piain FONLL-A-damp 



LHC 7 TeV - Different HQ sctiemes 



FQf^LL-A-piain FQt^LL-A-damp 



LHC14TeV - Different HQ si 
Zt^-VFN FQf^LL-A-piain FQt^LL-A-damp 



Figure 30: Graphical representation of the results of Table IHl 



5.5 Comparison with present and future HERA and data 

In conclusion, we look at NNPDF2.1 predictions for F2 and F^, which are especially 
sensitive to the treatment of heavy quark mass effects. For instance, heavy quark mass 
corrections to F^ for < 20 GeV^ are larger than 30%, see Fig. [3l In Fig. [36] we compare 
to the best fit result the HERA F2 data which have been included in the present analysis: 
in general the agreement is rather good, though the lowest and x bins have been 
removed from the fitted dataset, because O (a^) heavy quark corrections, not included 
in the present analysis, are large there [TT]. The NNPDF2.1 predictions for Fl{x,Q'^) is 
compared to pubhshed ZEUS [80j and HI [M] data in Fig. El Note that while the HI 
data are included in the fit, they have rather large uncertainties and thus carry very little 
weight in the global fit. Predictions obtained using NNPDF2.0 PDFs, but including heavy 
quark mass effects in the computation of the structure function through FONLL-A are 
also shown. This comparison is particularly interesting, because heavy quark mass effects 
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Figure 31: The parton luminosities Eqs. ((35ll36| for NNPDF2.1 compared to NNPDF2.0, CTIO 
and MSTW2008: from left to right and from top to bottom $cc, ^bb, ^bg- All 

Imninosities are plotted as ratios to the NNPDF2.1 central value. PDF sets with as{Mz) = 0.119 
have been used in all cases. All uncertainties shown are one-cr. 



LHC7TeVn LHCTTevH LHCyTevH 




Figure 32: Same as Fig. [311 but will all luminosities normalized to their respective central values. 



are quite large especially at low Q^: this correction is included here with both sets, though 
NNPDF2.0 PDFs were determined without it. The good agreement between results found 
using the two sets shows that NNPDF2.0 PDFs are quite accurate despite the lack of 
heavy quark mass corrections in the fit. 

We now turn to the predictions in view of upcoming combined HERA data. In Fig. [38] 
we provide the NNPDF2.1 predictions for in the kinematic region of the upcoming 
combined HERA data 0. The increase in uncertainty at small-x is driven by the larger 
uncertainty on the gluon at small-x, as seen in Fig. [T2j We also show the results using 

^We thank S. Glazov for providing us this information. 
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Ratio to NNPDF2.1 a,=0.1 1 9, Q' = 2 GeV^ | Ratio to NNPDF2.1 a,=0.119, tf = 2 GeV^ | 




Figure 33: Comparison between NNPDF2.1 sets with different values of the strong coupling, 
shown as ratios with respect the reference fit with = 0.119. The PDFs shown are the gluon at 
small and large- a; (upper plots), the triplet at large- 2: and the singlet at small- a; (lower plots). 

FONLL-A with both NNPDF2.0 and 2.1 input PDFs. The NNPDF2.1 results have been 
compared with preliminary combined HERA Fl dataset in Ref. [82]. 

In conclusion, we give predictions for F2 in Fig. [39] in the range of upcoming combined 
HERA datall These predictions are obtained with heavy quark mass effects included up to 
0(0^) through the FONLL-B scheme, but using input PDFs determined with 0{a) heavy 
quark mass corrections. The fact that FONLL-A and B coincide for moderate and large 
values of Q^, where the data included in NNPDF2.1 he, justifies the use of FONLL-B 
to extrapolate to the low region with the same input PDF set. PDF uncertainties at 
small-x and are rather large, suggesting that the combined HERA F2 data will impose 
severe constraints on the small-x gluon; comparison with preliminary data [S2] suggests 
that the NNPDF2.1 will be in very good agreement with the HERA data down to the 
smallest values of Q^. 



We thank K. Lipka for providing these plots. 
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NNPDF2.1 - Q^ = 2GeV^ NNPDF2.1 - Q^ = 10''GeV^ 




Figure 34: Correlation coefficient between PDFs and as{Mz) computed assuming as{Mz) = 
0.119 ± 0.0012 at 68% C.L. Results are shown at low scale (Q^ = 2 GeV^, left) and high scale 
(Q2 = 100 GeV^, right). 
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Figure 35: Determination of the Weinberg angle from the uncorrected NuTeV data [5S], with 
[S-] correction determined from NNPDF1.2, NNPDF2.0 and NNPDF2.1. The uncertainty shown 
on NNPDF2.0 and NNPDF2.1 is the one-a PDF uncertainty only. 
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Figure 36: Comparison between tlie HI and ZEUS F| data included in the present analysis and 
the NNPDF2.1 best fit. The data uncertainty includes statistical and systematic errors added in 
quadrature, while the theoretical uncertainty is the PDF uncertainty only. 



m 2010 measurement | 




Figure 37: Comparison between HI [SI] (left) and ZEUS [50] (right) Fl data and NNPDF2.1 
predictions. Predictions using NNPDF2.0 with the FONLL-A are also shown. 
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Kinematics of combined HERA measurement | 
0.35- 




Figure 38: The NNPDF2.1 predictions for Fl in the kinematics of upcoming combined HERA 
data. Predictions using NNPDF2.0 with the FONLL-A are also shown. 




Figure 39: NNPDF2.1 predictions for with heavy quark mass effects inchided up to 0{a^) 
through the FONLL-B scheme in the kinematic region of upcoming combined HERA data. 
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6 Heavy quark mass dependence 



We will now discuss the dependence of PDFs on the values of heavy quark masses rric 
and mi), and present some preliminary investigations on uncertainties related to heavy 
quark masses, along the lines of the recent detailed study by the MSTW group [83]; as 
an outcome of this analysis we will present NNPDF2.1 sets with varying rric and rrib 
masses. We will first discuss how different features of the NNPDF2.1 PDFs depend on 
the values of rric and m^. We then briefly investigate the dependence on rric of some LHC 
observables and summarize how uncertainties on heavy quark masses can be treated in 
the Monte Carlo approach. Finally, as an example of a possible application we evaluate 
the combined PDF+m^ uncertainty for MSSM Higgs boson production. 





mc [GcV] 


mb [GcV] 


NNPDF2.1 


1.414 


4.75 


NNPDF2.0 [10^ 


V2 


4.3 


cTio [rej 


1.30 


4.75 


MSTW2008 [S5^ 


1.40 


4.75 


ABKM09 [22, 


1.50 


4.50 


HERAPDFl.O ESI 


1.40 


4.75 



Table 10: The default values of the heavy quark masses used in NNPDF2.1 and in several recent 
PDF sets. 

6.1 Dependence of PDFs on heavy quark masses 

The default value of heavy quark masses used so far are summarized and compared to 
those of other PDF sets in Table [TOj The dependence of PDFs on the heavy quark masses 
is studied by repeating the NNPDF2.1 fit with different mass values. In particular, we 
have repeated the reference fit for charm quark masses of 1.5, 1.6 and 1.7 GeV as well as for 
bottom masses of 4.25, 4.5, 5.0 and 5.25 GeV. It is important to observe that at the order 
at which we are working, the perturbative definition of the heavy quark mass is immaterial: 
indeed different definitions (such as, for example, the pole and MS mass definitions) differ 
by terms of 0{as). However, we are including heavy quark mass corrections up to 0{as) 
only, so the difference is subleading (it becomes relevant once one includes 0{a1) heavy 
quark corrections, for example using the FONLL-B scheme). Therefore, the value of the 
quark mass in our PDF determination (as well as in other PDF determinations based 
on an NLO ACOT treatment of heavy quarks, such as CT/CTEQ) can be equivalently 
interpreted as, say, a pole mass or an MS mass. The MS mass is better known, and it has 
been recently shown ^1] to lead to perturbatively more stable results for deep-inelastic 
structure functions. 

Results are shown in Figs. I4U 14 II where the ratio of PDFs for different values of rUc and 
jTift to the reference NNPDF2.1 fit are plotted as a function of x for = 10^ GeV^. The 
dependence of the heavy quark PDFs on the value of the mass is easily understood: heavy 
quark PDFs are generated radiatively, and assumed to vanish at a scale equal to their 
mass. Therefore, a lower mass value corresponds to a longer evolution length and thus to 
a larger heavy quark PDF, and conversely. Thus, if one allowed [85] for an "intrinsic" [86] 
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Figure 40: Ratio of NNPDF2.1 PDFs obtained for different values of tlie cliarni quark mass to 
the reference NNPDF2.1 set at = 10" GeV^. Top left: gluon; top right: charm; bottom left: 
up; bottom right: down. 



heavy component (i.e. for a nonvanishing initial condition) this uncertainty would be 
absorbed in the initial intrinsic heavy PDF. Because of the momentum sum rule, if the 
charm PDF becomes larger, other PDFs are accordingly smaller (and conversely). For 
bottom in principle the same mechanism is at work, but in practice the effect on all other 
PDFs is negligible. 

6.2 Mass uncertainties and LHC observables 

The dependence of light quark distributions and the gluon on the charm mass displayed 
in Fig. Ho] is strong enough to affect the LHC standard candles studied in Sect. 15.11 at 
the percent level or more, as was discovered relatively recently [l5j. On the other hand, 
the dependence on the bottom mass of all PDFs but b itself is below the percent level, so 
only observables which depend on bottom are affected significantly: an example will be 
considered in Sect. [U^ below. 

Values of LHC standard candles computed as in Sect. 15.11 but using the PDF sets 
with different values of the charm mass of Fig. [JO] are collected in Table [TT\ and shown 
in Fig. I32J For completeness, we also give in Table [12] standard candles at 7 TeV for 
several values of the b mass. The variation of all standard candles is at the percent 
level for charm mass variations of order of 10%. It is interesting to observe that the 
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Q2 = 10" GeV^, ratio to NNPDF2.1 m = 4.75 GeV 



= 10" GeV^ ratio to NNPDF2.1 = 4.75 GeV 




Figure 41: Ratio of NNPDF2.1 PDFs obtained for different values of the bottom quark mass to 
the reference NNPDF2.1 set at = 10'' GeV^. Top left: gluon; top right: bottom; bottom left: 
up; bottom right: down. 



variation seen when modifying subleading charm mass terms is (recall Table [9]) of the same 
order of magnitude and in fact somewhat larger. This suggests that even though PDF 
uncertainties on standard candles are still dominant at present, theoretical uncertainties 
related to the treatment of charm will become relevant and possibly dominant as soon as 
PDF uncertainties are reduced by a factor of two or three. 

6.3 Combined PDF+m/i uncertainties and correlations 

Uncertainties which combine PDF uncertainty and heavy quark mass dependence are 
easily determined in a Monte Carlo approach, provided that PDF sets for several values 
of the quark masses are available. Here we provide several such sets, and more with a 
finer mass spacing and wider range will be made available in the future (such as already 
available for the MSTW08 sets [83]). Given sets of PDF replicas labelled by heavy quark 
mass values, such that PDF^'^'^'*'-'^ is the kij-th replica of the PDF set with heavy quark 

mass values rric^ and the mean value of any observable J- is 

(-^>rep = E E E -^(PDF(^'^-^),m«,mi^)) , (38) 

-''rep • 1 ■ 1 , -, ^ ^ 
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LHC 7 TeV 


W+Bi^ [nb] 


W-Bi, [nb] 


Z^B,j [nb] 


tt [pb] 


99 [pb] 


nic = 1.414 GeV 


5.99 ±0.14 


4.09 ± 0.09 


0.932 ± 0.020 


170 ±5 


11.64 ±0.17 


nic = 1.5 GeV 
nic = 1.6 GeV 
nic = 1.7 GeV 


6.06 ±0.17 
6.11 ±0.14 
6.14 ±0.14 


4.14 ±0.12 
4.17 ±0.10 
4.19 ±0.09 


0.943 ± 0.024 
0.951 ± 0.020 
0.956 ±0.019 


169 ±6 
167 ±6 
166 ±5 


11.65 ±0.25 

11.70 ±0.21 

11.71 ±0.22 


P [o-, mc] 


0.44 


0.41 


0.48 


-0.31 


0.16 



LHC 14 TeV 


W+Bi, [nb] 


W-Bi, [nb] 


Z^Bfi [nb] 


ti [pb ] 


gg^H [pb] 


mc = 1.414 GeV 


12.00 ±0.27 


8.84 ±0.17 


1.99 ±0.036 


946 ± 19 


37.50 ± 0.40 


mc = 1.5 GeV 
mc = 1.6 GeV 
mc = 1.7 GeV 


12.01 ±0.31 
12.24 ±0.28 
12.37 ±0.28 


8.94 ±0.22 
9.02 ±0.20 
9.10 ±0.18 


2.01 ±0.04 
2.03 ±0.04 
2.05 ±0.04 


942 ± 24 
939 ± 22 
935 ± 19 


37.62 ± 0.62 
37.90 ± 0.55 
38.15 ±0.58 


P [o-, mc] 


0.48 


0.50 


0.56 


-0.19 


0.41 



Table 11: LHC standard candles at ^/s = 7 TeV (upper table) and 14 TeV (lower table) obtained 
using NNPDF2.1 fits with different values of the charm mass rUc] the values in the top line of each 
table are the same given in Sect. 15. ll The bottom line of each table gives the correlation coefficient 
between the observable and the mass. 





W+Bi, [nb] 


W-Bi^ [nb] 


Z^Bu [nb] 


tt [pb ] 


99 [pb] 


nib = 4.25 GeV 
nib = 4.5 GeV 
nib = 4.75 GeV 
nib = 5.0 GeV 
mb = 5.25 GeV 


5.97 ±0.12 
5.95 ±0.21 
5.99 ±0.14 
5.99 ±0.12 

5.98 ±0.11 


4.07 ± 0.08 
4.07 ±0.11 

4.09 ± 0.09 
4.11 ±0.07 

4.10 ±0.07 


0.930 ±0.016 
0.928 ±0.025 
0.932 ± 0.020 
0.932 ±0.016 
0.930 ±0.015 


170 ±6 

171 ±7 
170 ±5 

170 ±5 

171 ±6 


11.58 ±0.26 
11.64 ±0.18 
11.64 ±0.17 
11.64 ±0.17 
11.66 ±0.18 



Table 12: LHC standard candles ai ^/s ^ 7 TeV obtained using NNPDF2.1 fits with different 
values of the bottom mass TOf,; the values in the top line of each table are the same given in 
Sect. O 



where 



N, 



rep 



i=i j=i 



(39) 



(i i) 

is the total number of replicas, and N^ep are distributed according to a two dimensional 
gaussian distribution with mean {m'^\m^^) and width {6mc,Snib) (assuming the values 
of charm and bottom masses are uncorrelated) : 



iV(«'i) 



rep' « exp 



V 



(i) (0) 
TUc — nic 

261 _ 



h 



(40) 



The combined PDF±m/j uncertainty is then the standard deviation of the observable over 
the replica sample: 

5pDF+m,-^ = J{T^)-{Tf , (41) 
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LHC7TBV 



LHC7TBV 



LHC 7 TeV 



4GeV I.SGeV 1 .6 GeV 



4GeV 1.5 GeV 1 .6 GeV 1 7 Ge*. 



4 GeV 1.5 GeV 1.6 GeV 1 7 Ge' 



NNPDF2.1 - LHC14 TeV - Char 




NNPDF2.1 -LHC14TeV-Ch: 




NNPDF2.1 - LHC14 TeV - Char 




Figure 42: Graphical representation of the results of Table [TTl 



where averages over replicas are to be understood as in Eq. ([38]) . Of course, a different 
probability distribution (possibly including a correlation between heavy quark mass values) 
could be assumed instead of Eq. (jlDj) . 

We can easily compute the correlation between PDFs and heavy flavour masses nih- 

(mhPDF (x, Q2) \ _ {rrih),,^ (PDF (x, Q^) ) 

f''mftO"pDF(x,Q2) 

where averages over replicas are to be understood in the sense of Eq. (|38p . The correlation 
Eq. (|i2]) . computed assuming rric = 1.55 ± 0.15 GeV and rrih = 4.75 ± 0.25 GeV, is 
displayed in Fig. [331 as a function of x for = 10^ GeV^. As discussed in Sect. 16. H 
as the mass is increased the corresponding heavy quark PDF is reduced, i.e. the PDF is 
strongly anticorrelated with its mass. As a consequence of the heavy quark suppression 
other PDFs are enhanced and this appears as a positive correlation, though the effect 
for bottom is negligibly small. Correlations between m/j and physical observables can be 
computed analogously. They are given in Tab. 1111 and can be used for a quick estimate of 
the corresponding uncertainty. 

6.4 The combined PDF+m;, uncertainties in MSSM bb ^ H production 

As an illustration of the procedure to combine PDF and nib uncertainties discussed in 
Section 16. 3|, we have evaluated the combined uncertainty for Higgs production in associ- 
ation with bottom quarks [87j. Higgs production via bottom fusion is enhanced in the 
MSSM in large tan /3 scenarios as compared to the SM, so this channel is important for 
supersymmetry searches. 

We have used the code of Ref. [87] to computed the bb ^ H cross-section to NLO in 
the MSSM, using the NNPDF2.1 sets with variable of Sect. 16.11 For other physical 
parameters we take the default values. Results are shown in Fig. [35] as a function of the 
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NNPDF2.1, Correlation PDFs and m^, = lo'* GeV^ 



NNPDF2.1, Correlation PDFs and m^, = lO'* GeV^ 




Figure 43: Correlation between PDFs and tire heavy quark masses at a typical LHC scale = 10^ 
GeV^: charm mass variations (left plot) and bottom mass variations (right plot). These correlations 
quantify the qualitative behaviour observed in Figs. HDl HTl 



Higgs mass for LHC 7 TeV, with two different uncertainty ranges for the bottom mass (the 
current PDG [6l] quotes ILq'qy as uncertainty on the MS b mass). Even with the smaller 
range the mass uncertainty is not negligible in comparison to the PDF uncertainty. 

It would be interesting to extend this analysis to several LHC processes which are 
expected to depend significantly on heavy quark masses, such as for instance t-channel 
single-top production. 




Figure 44: Combined PDF+mh uncertainties on the total cross section for bb H Higgs pro- 
duction at the LHC 7 TeV with — 4.75 GeV. The absolute cross section (left) and relative 
uncertainty (right) are shown. In the right plot, the PDF-only uncertainty is compared to the 
combined PDF+mi, uncertainty with Sm^ = 0.15 or Smt = 0-25 at one-cr. 
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LHAPDF name 


/vrmax 


mc [GeV] 


mt, [GeV] 


Qmatch [GeV] 


Ols (Qmatch) 




NNPDF21_100.LHgrid 


6 


V2 


4.75 






0.11900 


NNPDF21.FFN.NF3.100.LHgrid 


3 






V2 


0.359912 


0.10585 


NNPDF21.FFN.NF4.100.LHgrid 


4 


V2 




4.75 


0.218200 


0.11343 


NNPDF21.FFN.NF5.100.LHgrid 


5 


V2 


4.75 


175 


0.108283 


0.11900 



Table 13: NNPDF2.1 PDF sets with maximum fixed flavom' number of active quarks. In all cases, 
the number of flavours of the reference NNPDF2.1 PDF set is frozen at Qmatch' ^^^^ PDFs are 
then evolved upwards with a fixed number of fiavours. For each set, the values of the heavy quark 
masses, matching scale, and strong coupling at the matching scale and at = are shown. 



7 PDFs with Fixed Flavour Number 

The NNPDF2.1 PDFs discussed so far are determined in a factorization scheme in which 
the number of flavours depends on the scale, and in particular it varies from Nf = 3 
when < to Nf = 6 when > m^. Such a scheme is advantageous in that it 
includes terms to all orders in up to the desired logarithmic order (NLO in our case) 
both for light and for heavy quarks, while in a general mass scheme such as the FONLL- 
A scheme used here heavy quark mass terms are also included up to some fixed order 
in as- Indeed, given that the charm mass is of the same order as the starting scale of 
perturbative evolution at which PDFs are parametrized (in fact, they coincide for the 
default NNPDF2.1 set), the resummation to all orders in Og of charm mass logarithms of 
the form InQ^/m^ is as important as that for any other parton distribution. However, the 
LO and NLO resummation of logarithms related to heavier flavours is usually rather less 
important at scales relevant for LHC phenomenology, especially for top, but in practice 
sometimes also for bottom. In these cases, PDF sets in which the maximum number of 
flavors is fixed at some value lower than NJ^'^^ = 6 may lead to equally accurate results. 
Furthermore, use of these sets is necessary in conjunction with matrix elements computed 
with a number of active flavour smaller than six, such as single top production [88j, Higgs 
production in association with bottom quarks [89], as well as with Monte Carlo codes 
based on similar computations, such as the HVQDIS [90] Monte Carlo, widely used for 
the analysis of F2, which is based on the Nf = 3 scheme computation of the observable. 

With these motivations, we have constructed PDF sets in which the maximum number 
of flavours is Nf = 3, Nf = 4 and Nf = 5 (see Table [T3|) . These are simply constructed by 
freezing the number of flavours at some scale Qmatch which is thus viewed as a matching 
scale between a scheme in which the number of flavours depends on scale (for < Qmatch) 
and a scheme in which the number of flavours is fixed (for > Qmatch)- ^^^^ refer to 
these as Fixed-Flavour Number (FFN) PDFs, though this is strictly speaking a misnomer 
except in the Nf = 3 case: for instance, below the top threshold the default NNPDF2.1 
set and the Nf = 5 set are identical. Note that if these PDFs are to be used with 
matrix elements computed with the given number of fiavours, the strong coupling must 
be consistently determined with the same fixed number of flavour, lest a spurious Nf 
dependence be introduced in physical observables (see Ref. |9l], where the effect of an 
incorrect choice in this respect is also estimated). 

In order to illustrate the differences between these PDFs in Fig.HSjwe compare the FFN 
PDFs in the Nf = 3 and Nf = A to the reference NNPDF2.1 PDFs at the scale = 192 
GeV^. Differences stem both from the milder evolution due to the reduced number of 
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Ratio to NNPDF2.1,Q^ = 100 GeV' | Ratio to NNPDF2.1, = 100 GeV^ | 




Figure 45: Comparison between the default NNPDF2.1 set and tlie iV/ = 3 and Nf ^ A FFN 
scheme sets at the scale — 10^ GeV^. Results are shown as ratios to the default set. The singlet 
and the triplet (top) and the gluon in a linear and logarithmic scale (bottom) are shown. 

quark flavours and from the correspondingly smaller value of a^. This is particularly 
clear in the case of the singlet at small-x which is substantially smaller in the Nf = 3 
scheme due to the missing contribution from the charm and bottom PDFs. In the Nf = 4 
case differences are smaller both because we are closer to the heavy quark threshold and 
because now only bottom is not included into the beta function running and the DGLAP 
evolution equations. 

A similar comparison, but now as a function of scale for fixed x is performed in Fig. H6l 
The default and FFN PDFs coincide below the matching scale, and become increasingly 
different as grows. Differences are larger for the gluon, which is coupled by evolution 
to the singlet, which depends on A'^^ already at leading order. 

Several commonly used PDF sets, such as for instance MSTW08 [2T] and CTIO [T6] use 
Nf = 5 as a maximum number of flavours. In most cases this makes very little difference 
at LHC energies, but, again, care must be taken that a scheme in which Nf = 6 above top 
threshold is sometimes used matrix element calculations such as for example, the NLO 
Higgs production cross section of Ref. |92j. To illustrate the size of the effects involved, in 
Fig. [17]we compare the gluon and the singlet PDFs for the reference NNPDF2.1 set and 
the FFN set with Nf = 5 as a function of scale at x = 10~^. 

It is beyond the scope of this work to study the phenomenological implications of these 
FFNS PDFs. However several applications can be envisaged, such as helping to determine 
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Ratio to NNPDF2.1 GM-VFN, x=10-^ 



Ratio to NNPDF2.1 GM-VFN, x^lO"^ 



■b 



I NNPDF2.1 
NNPDF2.1 FFN t^=3 
NNPDF2.1 FFN Np=4 



Q' [ GeV^ ] 



Ratio to NNPDF2.1 GM-VFN, x=0.2 



1 



I NNPDF2.1 

NNPDF2.1 FFN Np=3 
. NNPDF2.1 FFN N.=4 



Q^t GeV^ ] 



10' 



10* 



I NNPDF2.1 
NNPDF2.1 FFN Np=3 
NNPDF2.1 FFN Np=4 



to' 

[ GeV' ] 



10= 



Ratio to NNPDF2.1 GiVI-VFN, x=0.2 



a 



















^1 NNPDF2.1 




NNPDF2.1 FFN Np=3 




NNPDF2.1 FFN Np=4 



•to' , 10' 



Figure 46: Comparison between the default NNPDF2.1 set and tlie iV/ = 3 and Nf ^ A FFN 
sets as a function of scale. Results are shown as ratios to the default set. The gluon and singlet at 
X = 0.01 (top) and the triplet and the total valence at a; = 0.2 (bottom) are shown. 



the optimal scales in the massless computations when comparing with the massive result, 
where in each case the corresponding PDF set in the same scheme has been used. We will 
discuss these applications in future work. 



Ratio to NNPDF2.1 GM-VFN, x^lO" 1 Ratio to NNPDF2.1 GM-VFN, x^lO " 1 




Figure 47: Comparison between the default NNPDF2.1 set and the Nf ~ 5 FFN set as a function 
of scale, shown as ratios to the default set. The gluon (left) and singlet (right) are shown. 
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8 Conclusions and outlook 



The NNPDF2.1 PDFs presented here fulfill most of the requirements of an ideal NLO 
parton set: they are based on a global dataset which includes most if not all relevant deep 
inelastic and hadronic data, they are free of parametrization bias, they are provided with 
reliable and statistically meaningful uncertainty estimates, they include NLO contributions 
without resorting to X-factors, they include a consistent treatment of heavy quark mass 
effects, and they are available for a variety of values of the strong coupling and heavy quark 
masses, which allows for the estimate of the associated uncertainties. We find that the 
heavy quark mass effects, while substantial for observables which probe directly the heavy 
quark distributions, are rather small for more inclusive observables. In particular the 
benchmark LHC cross-sections (inclusive vector boson, top and Higgs production) change 
by little more than one standard deviation when going from NNPDF2.0 to NNPDF2.1. 

An important missing feature of NNPDF2.1 is a reliable estimate of theoretical uncer- 
tainties related to missing higher order corrections (or equivalently, to choices of renor- 
malization and factorization scale). The obvious direction for future improvement is thus 
the determination of NNLO corrections to PDFs, the inclusion of mass efects at 0(q!^), 
and the inclusion of resummation corrections at large and small x. We expect all these 
corrections to be generally rather less than current PDF uncertainties: nonetheless, they 
need to be computed. 

To the extent that the criteria fulfilled by NNPDF2.1 are the dominant ones for accu- 
rate phenomenology, NNPDF2.1 is perhaps the most reliable parton set currently avail- 
able: no other set fulfills all of these criteria. We believe it to be adequate for precise 
phenomenology at the LHC, at least for the levels of experimental precision to be ex- 
pected in the near future. 



All the NNPDF2.1 PDF sets that have been discussed in this work are available from 
the NNPDF web site, 

http : //sophia . ecm . ub . es/nnpdf 

and will be also available through the LHAPDF interface [93j . 

The PDF sets that have been produced in the present analysis and that will be available 
in LHAPDF are the following: 

• The reference NNPDF2.1 sets, sets of N^-cp = 100 and 1000 rephcas: 
NNPDF21_100.LHgrid and NNPDF21_1000 . LHgrid. 

• NNPDF2.1 sets of A^'rep = 100 rephcas with as varied from 0.114 to 0.124 with steps 
of Sos = 0.001: 

NNPDF21_as_0114_100. LHgrid, . . ., NNPDF21_as_0124_100 . LHgrid 
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• NNPDF2.1 sets with combined PDF+ag uncertainty: 
NNPDF21_as_0119_pni_00012_100.LHgrid, 
NNPDF21_as_0119_pm_00020_100.LHgrid, 

NNPDF21_as_0119_pm_00012_50 . LHgrid and NNPDF21_as_0119_pm_00020_50 . LHgrid 
Note than in this case the corresponding NNPDF2.0 sets have also been produced. 

• NNPDF2.1 sets based on reduced datasets: DIS only, DIS+DY, DIS+jets: 
NNPDF21_dis_100 . LHgrid, NNPDF21_dis+dy_100 . LHgrid, 

and NNPDF21_dis+jet_100. LHgrid, as well as NNPDF21_dis_1000. LHgrid 

• NNPDF2.1 sets of A^rep = 100 replicas with varying charm mass: 
NNPDF21_mc_150_100 . LHgrid, NNPDF21_mc_160_100 . LHgrid, 
and NNPDF21jnc_170_100. LHgrid. 

• NNPDF2.1 sets of A^^cp = 100 replicas with varying bottom mass: 
NNPDF21_mb_425_100 . LHgrid, NNPDF21_nib_450_100 . LHgrid, 
NNPDF21_mb_500_100. LHgrid and NNPDF21_mb_525_100. LHgrid. 

• NNPDF2.1 sets in the Nf = 3, Nf = A and Nf = 5 FFN schemes: 
NNPDF21_FFN_NF3_100 . LHgrid, NNPDF21_FFMJJF4_100 . LHgrid, 
and NNPDF21_FFNJJF5_100. LHgrid. 
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A Mellin space implementation of neutral current structure 
functions 

In this Appendix we compute the analytic Melhn transform of the x-space O (os) heavy 
quark neutral current coefficient functions and we discuss the implementation and bench- 
marking of FONLL neutral current structure functions in the FastKernel framework. The 
corresponding results for charged current structure functions are collected in Appendix iBl 
The x-space gluon O {og) heavy quark coefficient function is given by Eq. Its 
Mellin transform is defined in the standard way as 



It is easy to see that the integral Eq. ()43p can be written in the following way 

1 + 1; 



V 



C7(y'^(Ar,e) = Tru^ dtt^~^ {[l + 2a{2e - l)t + 2a\l - Qe - 46^)1^] In ^ 



-[l + Aa{e-2)t-4a^(e-2)t^]v^ (44) 

= T^a^ rdtt^-i|[l + (l-3a)t-i(l + 4a-9a2)t2]inl±^ 
Jo ^ 1 — u 

- [l + (1 - 9a)t - a{l - 9a)t'^] , (45) 

where we have defined a (e) = (l + 4e)~^ to simplify the coefficients. The integrals we 
need are thus 

Ji{N) = r dtt^'Hn^-^ , J2{N)= [\tt^-\, (46) 
Jo 1 - ^ Jo 

since extra powers of t can be accommodated by a shift in N by an integer. Here as usual 

t; = (l-t)l/2/(l_ai)l/2^ 

The two integrals that we need are related by an integration by parts. To show this, 
we need 

d - 1 + V dv d , 1 + V 

— m = — — m 

dt 1 — V dt dv 1 — V 

1 1-a \f2{l-at) 



2(l-t)i/2(i_at)3/2A {l-a)t 
1 1 



Thus 



i (l-t)l/2(l_a^)3/2• 



^^ , 1 + 1 



(47) 



J.(iV) = -- d«"-,„_ = -/(Ar), (48) 



where we have defined 



-1 

/(iV)=/ dtt^-^{l-t)-^/^{l-at)-^/^. (49) 
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Note that the boundary term in the integration by parts vanishes for all ReA^ > 0, and 
thus its analytic continuation vanishes for all A'^, so it can be safely ignored. Trivially 
J2{N) = I{N) — I(N + 1). The integral I{N) may be evaluated in the usual way in terms 
of a standard hyper geometric function: 

rfiv)rfi) 

= r(Ar + |) 2FiiiN,N + l;a). (50) 

Note that when a = 0, this reduces to B{N, |), as it should. 
We thus get 

C(";)'^ {N, e) = T«a^{ [j^I(N) + ^/(iV + 1) - ll±jf=I^I{N + 2)] 
- [I{N) - I{N + !) + (!- 9a)(/(iV + 1) - /(iV + 2)) 
-a(l - 9a)(/(iV + 2) - I{N + 3))] } 

= r^a^{(i - l)/(iV) + + 9a)I{N + 1) 



' 1 l+4a-9a^ 
^2 N+2 



(1 + a)(l - 9a))I{N + 2) - a(l - 9a)I{N + 3)} . (51) 



This result is the required ingredient to implement the FONLL-A neutral current structure 
functions in the FastKernel framework. 

A cross-check of the Mellin transform of the massive coefficient function Eq. ()5ip is 
provided by the fact that its massless limit coincides with the Mellin transform of the 
x-space massive asymptotic {ni,0) coefficient function, Eq. ([8]). To this purposes, we need 
to expand Eq. (f5T|) near a = 1. Near a = 1, i.e. e = 0, we need to use the asymptotic 
expansion 

r(Ar + i) ~ r(n + i)r(Ar + n) 



[2V'(n + 1) - i^in + i) - V'(^ + n) - ln(l - a)](l - a)", (52) 



so that 



I(N) - 1 y. r(n+l)r(A^ + n) 

r(l)r(Ar) ^4^^ (n!)2 



[2V'(n + 1) - V(^ + ^) - i^i^ + n)- ln(l - a)](l - a)". (53) 

The n = term then gives the In e collinear divergence, which is subtracted by the massless 
coefficient function: as e — t- 

I{N) = - ln(4e) - 2je - H^) " ^i^) + 0{e). (54) 
Substituting in Eq. ([5T]) we obtain 

, r-A^3^3Ar2 - Ar(2 + Af + A^2wij^g^ ^^(jY))l 
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as expected: the coefficient of the singularity is precisely the LO anomalous dimension 
7gg^(A^). Therefore we have checked that the massless limit is properly reproduced, 

{N, e) + 0{e) = "O^'^ (iV, e) , (56) 

with the massive asymptotic coefficient function given by the Mellin transform of Eq. ([8]) , 
as desired. 

For completeness we also provide the correponding expressions for the O (a^) heavy 
quark coefficient function for the longitudinal structure function F^^c, which is implicitely 
contained in Eq. ([^ since F2^c = Ft,c + Fl,c- The x-space expression for the longitudinal 
heavy quark coefficient function is 



2. 1+^ 



-Sez log h ivz(l 

I — V 



(57) 



Its Mellin transform can be computed using the integrals discussed above, with the result 
ini),i f^r9l\- r^8ea ^^^^^^^ +4 {I{N + 1) - I{N + 2)(1 + a) + /(iV + 2)) 



The massless limits of the x- and A^-space results are straightforwardly computed and 
checked to be related by Mellin transformation as they ought to. 

Now we turn to the implementation and benchmarking of these results in the FastK- 
ernel framework. The major improvement in the FastKernel framework as compared to 
Ref. [TU] is the inclusion of heavy quark mass effects in deep-inelastic scattering structure 
functions, following the FONLL-A general-mass scheme [TT] . As discussed in Ref. (10) . 
FastKernel requires to write down all the DIS observables in Mellin space and precom- 
puting all the associated x-space Green's functions. Therefore, to extend FastKernel with 
FONLL structure functions we need to formulate FONLL in Mellin space. 

The x~space expression for the FONLL-A heavy quark structure functions, Eq. ([1]), 
can be easily written down in A^-space as follows: 

= F^\n,Q') (58) 

- e{Q'-m') {i-^)\f^-'\n,q')-f^^^'\n,q' 



with the default damping factor as threshold prescription. In order to implement Eq. ()58p 
in the FastKernel framework, we need the Mellin space expressions of the heavy quark 
coefficient function in the (n/), (n;, 0) and (n; + 1) schemes. While the last two are known, 
the former was not available in a closed form suitable for analytical continuation. The 
details of the computation have been presented above, and the desired result is Eq. dST]). 

With all the Mellin space heavy quark coefficient functions available, it becomes possi- 
ble to implement the FONLL-A heavy quark structure functions, Eq. (j58p into the FastK- 
ernel framework. To show that the A^-space implementation has the required accuracy, 
in Table for F2^c and in Table [T^ for F^^c, we compare the results for the Les Houches 
Heavy Quark benchmarks [53j for FONLL-A obtained with the FONLLdis code |lHI66j 
and with the FastKernel framework for various relevant values of Q^. The benchmark 
settings for the PDFs and are used for this comparison. What we can see is that 
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the accuracy on the FONLL heavy quark structure functions is essentially always below 
the percent level, enough for precision phenomenological studies. For completeness we 
also show the analogous results for the case of the massive scheme results, where similar 
accuracies are obtained. 





FONLL-A 


FFN 


X 


FONLLdis 


FastKcrnel 


Accuracy 


FONLLdis 


FastKernel 


Accuracy 


Q2 ^ 4 GcV^ 




0.1507 


0.1501 


0.4% 


0.1088 


0.1091 


0.3% 


10-4 


0.0936 


0.0931 


0.5% 


0.0697 


0.0698 


0.1% 


10-3 


0.0506 


0.0504 


0.4% 


0.0392 


0.391 


0.2% 


10-2 


0.0174 


0.0177 


1.5% 


0.0136 


0.0137 


0.7% 


Q2 = 10 GcV2 


10-^ 


0.563 


0.561 


0.4% 


0.3598 


0.3602 


0.1% 


10-4 


0.312 


0.311 


0.3% 


0.2007 


0.2011 


0.2% 


10-3 


0.1499 


0.1495 


0.3% 


0.0981 


0.0982 


0.1% 


10-2 


0.05056 


0.05052 


0.1% 


0.0328 


0.0327 


0.3% 


= 100 GeV^ 


10"^ 


2.28636 


2.28577 


0.02% 


1.9779 


1.9877 


0.5% 


10-4 


1.12186 


1.12082 


0.1% 


0.9161 


0.9184 


0.3% 


10-3 


0.48008 


0.47919 


0.2% 


0.3644 


0.3647 


0.1% 


10-2 


0.15207 


0.15200 


0.04% 


0.1037 


0.1038 


0.1% 



Table 14: Results of the benchmark comparison for the F2c{x,Q^) structure function in the 
FONLL-A scheme for the FONLLdis code [TT] and for the FastKernel code. Results are provided 
at the benchmark kinematical points in x, Q^. Results for the massive (FFN) scheme are also given 
for completeness. 





FONLL-A 


FFN 


X 


FONLLdis 


FastKcrnel 


Accuracy 


FONLLdis 


FastKernel 


Accuracy 


= 4 GeV2 


10' 


5 


0.0130174 


0.013094 


0.6% 


0.009077 


0.009081 


0.04% 


10- 


4 


0.008347 


0.008316 


0.4% 


0.005913 


0.005910 


0.05% 


10- 


3 


0.004795 


0.004778 


0.3% 


0.003511 


0.003509 


0.06% 


10- 


2 


0.001910 


0.001907 


0.2% 


0.001403 


0.001406 


0.2% 


02 = 10 GcV2 


10- 


b 


0.073235 


0.073022 


0.3% 


0.049856 


0.049982 


0.2% 


10- 


4 


0.041392 


0.041251 


0.3% 


0.028402 


0.028423 


0.07% 


10- 


3 


0.020754 


0.020707 


0.2% 


0.014463 


0.014456 


0.05% 


10- 


2 


0.007616 


0.007595 


0.3% 


0.005350 


0.005346 


0.07% 


g2 = 100 GeV2 


10- 


5 


0.471889 


0.4729 


0.2% 


0.3955 


0.397855 


0.6% 


10- 


4 


0.2236 


0.2235 


0.1% 


0.18656 


0.186914 


0.2% 


10- 


3 


0.0920 


0.09188 


0.1% 


0.0765 


0.076393 


0.1% 


10- 


2 


0.027822 


0.02782 


0.1% 


0.023079 


0.023100 


0.1% 



Table 15: Same as Table [TH for the Flc{x,Q^) structure function. 
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B Mellin space implementation of charged current struc- 
ture functions 



In this Appendix the analysis of Appendix A is repeated for charged current structure 
functions. 

The F2 charm structure functions in charged current DIS is given by Eq. ([9]). The 
expression for structure functions in neutrino-induced charged current scattering in the 
FFN scheme is 



2 27r 



+cf;)(e^/i^A)5(|„/i 



with i = 1,2,3. In Eq. (|59p have used the following definitions: 
s' = 2\V,s?s + 2|Krf|2[/d+ (1 - /)n]; / 



+ 



(59) 



(60) 



The explicit j;-space expressions of the O (a^) contributions C'i g(^) to the coefficient func- 
tions are given in Refs. |6HI62j. The standard structure functions are related to those 
defined in Eq. ([59]) through 



(61) 



so that 



Fl = F^ - 2xF{ = 2i [Tl - \Tl\ 



(62) 

Before Mellin- transforming the x-space quark coefficient functions of Refs. |61ll62j we 
rewrite them all in the form 



C^;;'^(x) = i^<5(l-x) + /(x) + [5(x)]+, 



where IT is a constant and f{x) is regular function in x G [0, 1] (so in general g{x) is not 
regular in x = 1). We get 



+ Cf 
+ Cf 



(1 + z^)lnz 



21 



\-z 



2(1 + z) ln(l - z) + (1 + z) ln(l - \z) + (3 - z) + 



z — z 



1 - \z 



In(l-z) ln(l-A2) 1 1 \-z 1 + , / 
4^ - 2—^^ - 2- + In A 



1 - z 



1 - z 



1-z 2 (1 - Az)2 1 - z 



(63) 
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+ Cf 



(-[ _|_ ^2-\ 1 / ? 

^ -2(l + z)ln(l-z) + (l + z)ln(l-Az)+ f2z + 2- - 



In(l-z) In(l-Az) 1 1 1-z 1 + z' 



1 - z 



+ 



l-z l-z 2(1-Az)2 l-z 



In A 



(64) 



+ Cf 



+ Cf 



2A 3 2A 

(1 + z^)lnz 



l-z 



2(1 + z) ln(l -z) + {l + z) ln(l - Az) + (1 + z) + 



In(l-z) In(l-Az) 1 I l-z l + z 



1 - z 



+ :t 



l-z l-z 2(1-Az)2 l-z 



In A 



with Ka = {1- A) ln(l - A)/A. 

The gluon coefficient functions do not need any further work and are given by 

C^^'\z) = rj(2z2-2z + l){21n(l-z)-21nz-ln[A(l-A)]} + 



(65) 



[4-4(l-A)]z(l-z) + (l-A) 
2(1 -A) 



+ 



zin / -2Az2ln 



1 - Az 
1 - Az 



(1 - A)z 



(1 - A)z 



1; 



(66) 



Tf{2z^ - 2z + 1) {2 ln(l - z) - 2 In z - ln[A(l - A)]} + 



[8 - 18(1 - A) + 12(1 - A)2]z(l - z) + (1 - A) 
6A(1 - A) 



1 



1 — Az , o , 1 — Az 
z In 2Az^ In ■ 



(l-A)z 



(l-A)z 



1 - Az 
1; 



+ 



C''^^\z) = r/(2z2-2z + l)<{21n(l-z)-21n(l-Az)+ln 



I: 



2(1 - A)z(l - z) + 2(1 - A) 



+ 



(1 + A)z2 In 



z In ■ 



1 - Az 



1 - Az 



(l-A)z (l-A)z 



(67) 



(68) 



In order to transform to the A^-space the above x-space expressions, in Tables [T6l and [T71 
we tabulate the Mellin transforms of all terms involved. In these tables we use the analytic 
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M[/](iV) 


6(1 -z) 


1 


z' 


1 

N + 1 


{l + z^)lnz 
l-z 




(l + z)ln(l-z) 


Si Si 1 
N N + 1 (A^ + l)2 


(1 + z) ln(l - Xz) 


2i^i(l,A^ + l,A^ + 2;A) lii(l-A) 

N{N + 1) N 

2Fi{l,N + 2,N + 3;X) ln(l - A) 
(A^ + l)(iV + 2) N + 1 




z - z'^ 




2Fi{l,N + 1,N + 2,X) 2Fi{l,N + 2,N + 3,X) 


1-Xz 


N + 1 N + 2 








^2Fi{l,N-l,N,X) 2Fi{l,N,N + l,X) 


1-Xz 


N-1 N 

2Fi{l,N + l,N + 2,X) 
N + 1 




1 - z 




2Fi{l,N,N + 1,X) 2Fi{l,N + 1,N + 2,X) 


1 - Xz 


N N + 1 




"ln(l - z)' 

l-z 


+ 






ln(l - Xz) 
l-z 


+ 


•^a(^) = Et Si{N + k)-Si{k)-^ 

k=i 




1 

l-z_ 


+ 






l-z 


+ 


2Fi{2,N,N + 2,X) A + ln(l-A) 


.(1- 


N{N + 1) A2 




"1 + ^2" 

_l-z _ 


+ 


AT N + 1 ^^'^2 



Table 16: Mellin transforms of the terms involved in the NLO charged current quark 
coefficient functions. 



68 



M[f]iN) 



2 In 



[z2 + (l-z)2]x 
1 - Z 



4 - 2N{N - 3) - N{N^ + N + 2){25i + ln[A(l - A)]} 
iV2(iv + i)(iV + 2) 



ln[A(l - A)] 



[z2 + (l-z)2)]x 



21n( ^— ^ ) +ln ' 



2 /A 



2A 



+ 



2Fi(l, N + 1, N + 2, X) 
A N + 1 ' N + 2J N + 1 

4(A-1) {N^ + N + 2){2Si-ln{^)] 



X{N + l){N + 2) 



N{N + 1){N + 2) 



zln 



I- Xz 



2Fi(l,iV + l,jV + 2;A) 
(A^ + l)2 



In 



I- Xz 

O^xfz 



2Fi{l,N + l,N + 2;X)-l 
X{N + l){N + 2) 



1 



l-Xz 



2Fi{l,N,N + 1;X) 
N 



1-Xz 



2Fi{l, N + 1, N + 2; X) 
N+1 



Table 17: Mellin transforms of the terms involved in the NLO charged current gluon 
coefficient functions. 



continuation of the harmonic sum 

fc=i ^ 

where ({I) is the Riemann ^-function, with ^(1) = ^em, fpiil — 1)) is the polygamma, and 
2Fi{a,b,c; N) is the Gauss hypergeometric function . 

As an example of use of Tables [TU1I171 we present here the complete A^-space quark 
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and gluon coefficient functions for 

^2 



Cf 



Si 



+ 



1 



1 /Si 
(Af + 1)2 + ViV^iV + 1 ' (iV + 1^2 



+ A 



2Fi(l,iV + l,iV + 2;A) 



iV(iV + l) 



In(l-A) , . 2Fi(l,iV + 2,Ar + 3;A) , ln(l - A) , 



N 



+ 2 



(iV + l)(iV + 2) ' iV + 1 
2Fi(l,iV- l,iV,A) 2Fiil,N,N + 1,X) 



2 2 



N -1 N -I 

2Fi{l, N + 1, N + 2, X) 
N + 1 

I (2Fi{2,N,N + 2,X) A + ln(l-A) 



+ 2[Sl + S2 



N 



2JxiN) - 2 



N 



N{N + 1) 



+ 



A2 



(69) 



4 - 2N{N - 3) - iV(iV2 ^ ^ ^ 2){25i + ln[A(l - A)]} 
' 2Af2(7v + l)(Ar + 2) 



+ 



8- 18(1 -A) + 12(1 -A)^ (1- A)2Fi(l,iV,iV + l;A) 
{N + l){N + 2) N 

2Fi{l,N +1,N + 2;X) 22Fi(l, TV + 1, iV + 2; A) - 1 



6A(1 - A) 



(iV + 1) 



(A^ + l)(iV + 2) 



(70) 
1 

N' 



As a cross-check of the Mehin space results, it is possible to compute the asymptotic 
limit A —7- 1 of these expressions. We need the asymptotic expansion of the hyper geometric 
functions, Eq. (j52p up to O (A — 1) terms. In particular, 

2Fi{l,N + 1,N + 2;X) = -(1 + N) (in {I - X) + je + + l)) + O ((A - 1)) , (71) 



2Fi(2, N,N + 2;X) = -N{1 + N) [in (1 - X) + + V-^"^ (N)) + O ((A - 1)) . (72) 

Substituting in Eq. (j69p . one can see that all collinear heavy quark logarithms and that 
the massless limit of the massive charged current heavy quark coefficient functions reduces 
to the usual ZM-VFN result, as we know from x-space. 

Now we turn to discuss the implementation and benchmarking of the above results 
into the FastKernel framework. Analogously to the neutral current sector, the FONLL-A 
charged current structure functions in Mellin space can be written as 



^CC,FONLL^^^^2) 



(73) 



m 



2\ 2 
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with i = 1,2,3. The Melhn space expressions of the massive heavy quark coefficient 
functions have been computed above, and the other ingredients of Eq. (j73|) are their 
massless limits and the standard Melhn transform of the ZM-VFN coefficient functions. 

With these results, we have implemented the FONLL-A charged current structure func- 
tions Eq. ()73p into the FastKernel framework. As it has been done in the case of neutral 
current observables, here we benchmark the accuracy of this FONLL scheme implementa- 
tion. We use again the same settings of the Les Houches heavy quarks benchmark study. 
The benchmarking of the FONLL-A CC structure function implementation in FastKernel 
is performed for the charm production cross section in neutrino induced DIS, defined by 
Eq. (|34p ■ that combines all three charged current structure functions. We have checked 
that the comparison of individual structure functions has a similar level of accuracy. 

Results for the benchmark comparison are shown in Table [THl As discussed above, the 
FONLL-A calculation of charged current structure functions has been implemented in a 
x-space code, FONLLdisCC, that we will use for the benchmarking with the FastKernel 
implementation. Results are shown for various values of relevant for the analysis 
of experimental data. The accuracy is similar to the one achieved for neutral current 
structure functions (see Tables [TilllSp . at the per mil level, suitable for precision PDF 
determinations. 





FONLL-A 


FFN 


X 


FONLLdisCC 


FastKernel 


Accuracy 


FONLLdisCC 


FastKernel 


Accuracy 


= 4 GeV^ 


10-'^ 


163.14 


164.06 


0.6% 




158.70 


158.15 


0.3% 


10--* 


109.48 


109.55 


0.1% 




106.81 


106.64 


0.2% 


10-3 


69.24 


69.35 


0.2% 




67.86 


67.88 


0.1% 


10-2 


37.75 


37.87 


0.3% 




37.27 


37.30 


0.1% 


10-1 


13.56 


13.57 


0.1% 




13.53 


13.51 


0.1% 


= 10 GeV^ 


10-'^ 


279.31 


278.71 


0.2% 




261.49 


261.55 


0.02% 


10--* 


167.02 


166.85 


0.1% 




157.27 


157.11 


0.1% 


10-3 


92.90 


92.87 


0.03% 




88.33 


88.12 


0.2% 


10-2 


44.92 


44.93 


0.02% 




43.36 


43.23 


0.3% 


10-1 


14.50 


14.48 


0.1% 




14.26 


14.28 


0.1% 


Q2 = 100 GcV2 


10-'^ 


674.55 


674.53 


0.02% 




651.21 


645.94 


0.1% 


10-4 


345.73 


345.81 


0.02% 




331.17 


329.14 


0.5% 


10-3 


161.70 


161.78 


0.05% 




153.94 


152.36 


0.1% 


10-2 


64.20 


64.26 


0.1% 




61.11 


61.06 


0.1% 


10-1 


15.79 


15.83 


0.2% 




15.33 


15.42 


0.1% 



Table 18: Results of the benchmark comparison for the dimuon charm production cross section 
Eq. (p4)) . in the FONLL-A scheme for the FONLLdisCC charged current code and for the FastK- 
ernel framework. Results are provided at the benchmark kinematical points in x,Q^. Results for 
the massive (FFN) scheme are also given for completeness. The inelasticity variable in the dimuon 
cross section for this benchmark table has been taken to be y — 0.5. The Les Houches Heavy 
Quark benchmark settings (53j have been used for the comparison. 
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